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Abstract 

In this work, we investigate the phenomenology of minimal four-family MSSM 
supergravity theories containing an additional generation (t', b', t', u') of heavy 
fermions along with their superpartners. We constrain the models by demand- 
ing: gauge coupling constant unification at high energy scales; perturbative values 
for all Yukawa couplings for energy scales up to the grand-unification scale; ra- 
diative electroweak (EW) symmetry breaking via renormalization group evolution 
down from the grand-unification scale; a neutral LSP; and consistency with con- 
straints from direct searches for new particles and precision electroweak data. The 
perturbative constraints imply a rather light fourth-family quark and lepton spec- 
trum, and tan/? ^ 3. The lightest CP-even Higgs mass receives fourth-family loop 
corrections that can result in as much as a 30% increase over the corresponding 
three-family mass value. Significant fourth-family Yukawa coupling contributions 
to the evolution of scalar masses lead to unexpected mass hierarchies among the 
sparticles. For example, the r[ is generally the lightest slepton and the lightest 

squark is the b^. A significant lower bound is placed on the gluino mass by the 
simple requirement that the r[ not be the LSP. Sleptons of the first two families 
are much more massive compared to the LSP and other neutralinos and charginos 
than in the three-family models; in particular, all sleptons belonging to the first 
three families could easily lie beyond the reach of a \/s = 500 GeV e + e~~ collider. 
Consistency tests of the RGE equations via mass sum rules and relations are ex- 
plored. Relations between slepton masses and gaugino masses are shown to be 
very sensitive to the presence of a fourth generation. The most important near- 
future experimental probes of the four-family models are reviewed. A scenario 
with mt ~ mw and t — > £ix5 is shown to be inconsistent with universal soft- 
SUSY-breaking boundary conditions. Full four-family evolution of a s is shown to 
lead to a significant enhancement in inclusive jet and di-jet spectra at Tevatron 
energies when all sparticle masses are near their lower bounds. 



1. Introduction 

Despite the success of the Standard Model (SM), it is almost certainly in- 
complete. In particular, the full theory should include a quantum theory of 
gravity. However, the huge discrepancy between the characteristic Planck scale 
(Mpi ~ 10 18 GeV) and the weak scale, 0(mz), is not easily bridged without en- 
countering problems related to fine-tuning and gauge-hierarchy. Supersymmetry is 
currently the only fully quantitative and consistent resolution to these problems. In 



particular, the Minimal Supersymmetric extension of the Standard Model (MSSM) 
is an extremely attractive contender for physics beyond the SM due to its natu- 
ral prediction that the gauge couplings all unify at a common grand-unification 
(GUT) scale, Mjj. The precise particle content of the MSSM (comprising sparticle 
partners for all SM particles and exactly two Higgs doublet fields Hi, H2 and their 
higgsino partners) is crucial to this success. In addition, the MSSM provides an at- 
tractive explanation for the long proton decay lifetime and an attractive candidate 
(the lightest supersymmetric particle, or LSP) for cold dark matter, easily accom- 
modates the absence of significant flavor- changing neutral currents (FCNC's), and 
is completely consistent with current precision electroweak data. Of course, exper- 
iment has yet to reveal any evidence for SUSY, in either the sparticle or the Higgs 
sector. 

Unfortunately, supersymmetric models in general, and the MSSM in particu- 
lar, provide no explanation for the number of families (we denote the number of 
families/generations by N g ) or the patterns of fermion masses. Certainly, a super- 
symmetric extension of the standard three-family model is entirely consistent, but 
the possibility of incorporating one or more additional generations into the MSSM 

framework deserves exploration. In a previous work) 11 we investigated the minimal 
four-family, gauge-unified MSSM, with a large Dirac mass for the so as to be 
consistent with LEP constraints.' 31 

In particular, we explored the constraints upon the model obtained by gener- 
alizing the usual perturbative limit on the top quark Yukawa coupling to the 
requirement that all Yukawa couplings remain perturbative at energy scales below 
M\j . Assuming that the CDF and DO 'top-quark' events arise from a third genera- 

tion top with m t ^ 155 GeV, we demonstrated that perturbative Yukawa behavior 
requires a rather light fourth-family spectrum, namely m t i y ^ m t . The fourth- 
family lepton sector is even more strongly constrained: m T iy ^ 85 GeV, the upper 
limit occurring when int^rny are just beyond the reach of LEP. (That is, the upper 
limits for the fourth-family leptons versus quarks have a strong inverse correlation.) 
Additional results included: i) limits requiring tan/3 (tan/3 = i^/fi, v\ and t>2 be- 
ing the vacuum expectation values of the H\ and H2 scalar fields) to have a very 
modest value, roughly 1 ^ tan/3 ^ 3, in order to avoid perturbativity problems 
for the Yukawas; ii) the impossibility of imposing the X^Mjj) / \ t {Mjj) = 1 bound- 
ary condition in four-family models; and iii) an increase in the predicted value for 
as(mz) by ~ +3% from the purely gauge coupling contributions to gauge cou- 
pling running, which could, however, be largely compensated by two-loop Yukawa 

contributions to gauge coupling running.' 11 

In this paper we elaborate upon our earlier results and extend our four family 
study to include the superpartners of the fourth-family quarks and leptons. We 
adopt the conventional framework in which the MSSM parameters are determined 
in the context of minimal supergravity by universal (at Mjj) soft-SUSY-breaking 

* Wc define the fourth generation by the CKM matrix hierarchy \V t b\ 2 , \ Vvv | 2 ~ b with 

|^ 6 | 2 ,|vw| 2 «i. 
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masses for the gauginos (mjrt) and for the squark and Higgs boson fields (mo), and 
universal soft Yukawa couplings (A); of course, these typically evolve to (rather 
disparate) weak-scale values. Additional crucial parameters are: tan/9 (defined 

above) and /i, the coefficient of the H1H2 mixing term in the superpotential. Aside 
from having considerable theoretical motivation, this approach has the advantage 
that relatively few independent parameters, tan f3, mt, mo, A, (along with the 

sign of //), are sufficient to completely specify the theory at the weak scale.^ 

Specific supergravity (SUGRA) and string models make definite predictions 
for the relative sizes of mo, A and m^ 2 - It is convenient to specify a given model 

in terms of the ratios £4 = A/mi/2 and £0 = fflo/%/2- Two models provide 
particularly useful benchmarks. The first is the string-motivated dilaton boundary 

condition scenario specified by £0 — l/\/3, £a — ~ 1- This boundary condition set 
emerges universally in all string models where supersymmetry breaking is dilaton 
dominated. It represents a middle-of-the-road choice in that the gaugino mass m 1 / 2 
and the soft scalar mass mo are both of importance in the final low-energy values of 
the squark and slepton masses, but sleptons are generally significantly lighter than 
squarks, and both are generally lighter than the gluino. A more extreme boundary 
condition choice is the 'no-scale' model with £0 — £,A — 0, which can also arise in 
certain string and supergravity approaches. In this model, supersymmetry breaking 
arises entirely from the gaugino mass m 1 / 2 at Mjj, with all other supersymmetry 

breaking parameters generated by RGE evolution as the energy scale decreases. 
Sleptons are still lighter compared to squarks, and the gluino mass is generally the 
largest. A brief review of these two boundary conditions can be found in Ref. [5]. 
As a final benchmark possibility we shall also consider £0 = —£4 = 1. In this case, 
the slepton and squark masses turn out to be large (often larger than the gluino 
mass) and fairly similar in size due to the dominance of the mo source term. 

Our study is designed to complement the existing MSSM and supergravity 
(SUGRA) studies, virtually all of which have assumed three generations) 61 and the 
many earlier studies of a four-family SM. Since the LEP experiments rule out 

rgi 

the possibility of an additional new sequential 'light' neutrino, the fourth-family 
neutrino must be quite massive, m v i > 45 GeV, and the fourth family would seem 
to not be truly 'sequential'. However, there are many indications of small non-zero 
mass for the neutrinos of the first three families, in which case the much larger mass 
for the v' is no different than the large value of the top quark mass as compared 
to the masses of the other quarks. Both of these large generational hierarchies 
must find explanation in physics beyond the MSSM. A fourth family with a heavy 

neutrino is no more unnatural than a third family with a heavy top quark.' 91 

In the minimal supersymmetric model with N g < 4, the gauge couplings 

unify perturbatively at a common scale M v ~ 2 - 5 x 10 16 GeV. However, 
as noted earlier, for N g — 4we must relax the often- imposed theoretical prejudice 

f For a review of this approach, sec Ref. [4] and references therein. 
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that \b{Mu) = \ t (Mjj) if we are to generate an acceptable prediction for mfc/m r 
at low energies. We do not regard this as a significant difficulty since there are 
many reasons why the Yukawa couplings might not be unified at Mjj. In addition 

to potentially large weak and GUT-scale threshold effects) 111 a different field con- 
tent can drastically modify the GUT-scale relations among the Yukawa couplings. 
For example, if a 45-Higgs is added to the SU(5) GUT theory, then the strict 

relation \(Mjj) / 'X r (Mjj) = 1 is no longer valid.' 121 More recently, the authors in 
Ref. [13] have shown that if the theory contains additional heavy fermions at Mjj 
then \b{Mu) / \ T {Mu) < 1 can naturally result. 

Thus, it is of considerable interest to study the phenomenology of a four-family 
MSSM model. In our study, we will delineate how current experiments can either 
eliminate or confirm the existence of a fourth family. We shall also discuss a num- 
ber of theoretical subtleties that arise in the implementation of four generations in 
the context of the SUGRA framework. Aside from updating the constraints on mf, 
my, and tan/3, we will particularly address the following specific issues, i) When 
all basic constraints are applied, what are the experimentally allowed regions of pa- 
rameter space? ii) Does the electroweak (EW) radiative breaking mechanism work 
in the N g = 4 case? iii) What are the additional radiative corrections to the light- 
est SUSY Higgs mass? iv) How is the general spectrum of the sparticles affected 
by the presence of four families? v) In particular, are there new constraints on 
MsuSY arising from fourth-family sparticle mass constraints? vi) More generally, 
how does the fourth family sparticle spectrum compare to the spectra of the first 
three families? And, vii) what experimental constraints are imposed by the latest 
DO and CDF top quark searches, as well as a global fit to the latest EW precision 
data? We examine each of these issues, although not precisely in the above order. 

Section 2 is devoted to a study of gauge and Yukawa coupling unification. In 
particular, we refine our earlier analysis to include exact (numerical) solution of the 
coupled two-loop gauge and Yukawa couplings. The parameter space regions given 
by demanding perturbativity for the Yukawa couplings and consistency with non- 
observation of sparticles and fourth-family fermions at LEP are specified. Section 
3 discusses the radiative electroweak symmetry breaking mechanism. Section 4 
describes the physical Higgs boson masses and their phenomenology in the four- 
generation MSSM. Section 5 discusses a constraint on the SUSY sparticle mass 
scales that is peculiar to the four-generation MSSM model. Section 6 presents, for 
a number of typical models, the sparticle mass spectra that arise from a four-family 
scenario, and delineates allowed regions of soft-supersymmetry-breaking parameter 
space after all direct experimental constraints have been imposed. Key model- 
independent features of the sparticle spectra are identified. Section 7 reviews the 
latest direct collider limits on a fourth family, including the implications of the 

latest DO and CDF data. There, we demonstrate that a scenario' 141 in which the 

top is light, mt ~ mw, but not observed because it decays by the mode t — > tix? 
(the CDF/DO events coming from t' — > bW in this scenario), is not consistent 
with universal soft-SUSY-breaking boundary conditions by virtue of the Section 5 

constraint which forces the t\ to be very heavy. In addition, a global fit to the latest 
precision LEP data is presented and the ensuing constraints on a fourth family are 
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discussed. We also explore implications of a fourth family for the inclusive-jet and 
di-jet spectra of the light quarks following from the dramatic slow-down in the 
evolution of a s once energy scales above Msusy are reached. Section 8 presents 
our summary and concluding remarks. The renormalization group formulae and 
beta functions that are frequently referenced in the text appear in Appendix A. The 
question of the accuracy with which tan/3 can be determined by squark/slepton 
mass measurements is discussed in Appendix B. 



2. Gauge and Yukawa Coupling Unification for N g = 4 

Gauge coupling unification and the low energy prediction for m&/m r (assuming 
\b(Mjj) = X t (Mjj)) at the two-loop level in an MSSM four-family model were 
first considered in Ref. [15]. For complete family representations, gauge coupling 
unification implies values for a^{mz) and Mjj that are independent of the number 
of families at the one-loop level. At two-loops, there is a weak dependence on 
N g . Fig. 1 shows the unification of the couplings in both the three- and four- 
family scenarios including two-loop gauge contributions to the the beta functions 
but not including the two-loop Yukawa contributions. Assuming the published 

value of a~m(mz) = 127.9 (in the MS scheme), and adopting sin 2 8w(mz) = 
0.2316 (see the discussion of the next paragraph), and a single SUSY breaking scale 
Msusy = m z * we predict a%{m z ) = (0.1283,0.1326), afj = (0.0432,0.0917) and 
Mjj = (2.98,5.73) x 10 16 GeV for N g = (3,4), respectively. Here, the superscript 
'0' indicates that the two-loop Yukawa effects have not been included. We see that 
there is a slight shift in 0:3 (mz), a significant shift in Mu, and a factor of two 

increase in a v as one moves from N g = 3 to N g = 4. For the more recent value 
of a~^{mz) = 129.08 (= 128.05 in the MS scheme), and adjusting sin 2 9w( m z) 
downwards by —0.0002, to 

sin 2 Mm z ) 

= 0.2314, the resulting predictions are 
al(m z ) = (0.129,0.133), cfy = (0.0432,0.092) and M v = (3.25,6.28) x 10 16 GeV 
for N g = (3, 4). In either case, inclusion of a fourth family raises a^{mz) by about 
3% when calculated at two-loop order without including Yukawa contributions to 
the two-loop beta functions. We shall shortly return to this issue. 

In all the calculations that follow, we shall employ the published value for the 
electromagnetic coupling of a~m{ m z) — 127.9 in the MS scheme and keep the MS 
value of sin 2 Bw^ m z) fixed at sin 2 9w( m z) = 0.2316. Holding sin 2 9w( m z) fixed is 
an approximation. The actual 'best-fit' value of sin 2 0\y(mz) for a given set of ex- 
perimental data depends upon the top-quark mass, the masses of the fourth-family 
fermions, the mass of the light SM-like Higgs boson of the model, and the masses 
of the superpartners of all four families (see, for example, Refs. [16,18]). If exact 
coupling constant unification is demanded, and definite boundary conditions for 
the soft-supersymmetry-breaking parameters are specified at the unification scale, 



* As described in Ref. [16], the effective Msusy is generally not far from mz even when some 
superpartners are heavy. 
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the appropriate 'best-fit' value for sin 2 6w( m z) value within the four-generation 
model could be obtained by a self-consistent iterative procedure such as described 
in Refs. [16,18]. We have chosen to avoid this complexity for this first study of 
the four-generation supersymmetric model. We also note that we ignore the dif- 
ferences between DR and MS couplings and masses. In particular, although it 
is the DR couplings that are most naturally required to unify in supersymmetric 
models, the differences between DR and MS couplings are not significant at scale 
Msusy ~ m z compared to other uncertainties. As described later, we allow a 
certain level of 'error' in the unification of the coupling constants in order to ac- 
count for the remaining experimental uncertainty in sin 2 Ow( m z), the slight DR 
vs. MS mismatch and the (small) variation in sin 2 6\y(mz) that would occur as 
we change 4th-family masses and SUSY parameters. We will see that the allowed 
parameter space is only slightly sensitive to a relaxation of exact unification, and 
certainly the general phenomenological features and issues that we discuss would 
not be significantly altered by a more precise treatment. 

The my , my parameter space for ttv = m T > = 50 GeV (consistent with LEP 

limits of 45GeV M ) and mt = 165 GeV, equivalent to mtipole) = 175 GeV, that 
is allowed after demanding perturbativity for the Yukawa couplings for all energy 
scales up to Mjj, as well as coupling constant unification at Mjj, including two-loop 
Yukawa contributions to the gauge coupling evolution equations, is illustrated for 
tan/3 = 1.5 and tan/5 = 2.2 in Fig. 2. Our precise criterion for perturbativity is 
that all Yukawa couplings obey < 3.3; this value ensures that one-loop effects 
dominate over two-loop effects!" 1 Two possible levels of gauge unification are con- 
sidered: unification of as with a\ and «2 to within 2.5% and to within 0.01%. 
Fig. 2 shows the allowed parameter space regions for these two cases. We observe 
that some of the excluded points on the border become allowed if the precision 
demanded for unification of the couplings is relaxed by even a few percent. Our 
procedure and these results will be explained in more detail shortly. 

The allowed parameter space regions illustrated in Fig. 2 are even more re- 
stricted than those given in our earlier work, Ref. [1], at these same tan/5 values, 
for two reasons: i) the higher value of m t (m t ) = 165 GeV (vs. m t (m t ) = 160 GeV) 
has been chosen so as to yield mtipole) ~ 175 GeV, consistent with the latest CDF 
and DO experimental results' 201 (the precise pole mass value depends on the spar- 
ticle spectrum, but only weakly); ii) the two-loop Yukawa effects in the running of 
the gauge couplings feed back into the Yukawa couplings themselves so that they 
violate perturbativity more easily. This latter point can be understood by noting 
that inclusion of Yukawa couplings reduces the value of 0:3 (Q) when mt^m^ have 
values near the boundary of the perturbatively allowed region. The reduction in 
0:3 (Q) in turn slightly reduces the magnitude of the negative gauge contribution to 
the one-loop component of the Yukawa beta functions as Q runs from mz to Mjj 
(see Eqs. (A.l), (A. 5) and (A. 6) in Appendix A), constraining the weak-scale value 



f Unless explicitly indicated, up until Section 7 all masses are the running m(m) masses, and 
not the pole masses. 
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of the Yukawas to be somewhat smaller than before in order to remain perturbative 
at Mjj. 

This reduction in a% for points at the parameter space edge is illustrated in 
Fig. 3 which shows how much the values of Mjj and as(mz) can be affected by 
including the large Yukawa contributions to the two-loop gauge coupling beta func- 
tions in the four-family scenario. The values of Mjj and a^{mz) are presented as 

contours in the the mf,m^ plane. The shifts in as(mz), Mjj from a®(mz) : Mjj 
agree well with our approximate solutions from Ref. [1]. For mt^m^ near their 
maximal allowed values, a^mz), Mjj can be lowered by as much as 5%, 16% re- 
spectively. We see that the N g = 4 values of OL^{mz) are essentially the same as 
the Ng = 3 value when large Yukawa effects are included at two-loop. 

The precise boundary of the perturbatively allowed region depends upon the 
precision demanded for coupling constant unification. In Fig. 2 we showed the 
allowed regions obtained by demanding unification of the coupling constants to 
within 2.5% or 0.01% (the latter being essentially equivalent to exact unification). 
The allowed region in the 0.01% case is slightly reduced compared to the 2.5% 
case. This sensitivity is implicitly present even before the two-loop Yukawa terms 
are included in gauge coupling evolution. As noted above, by increasing 0:3 the 
Yukawa blow-up is delayed. So if a^{Mjj) is allowed to be slightly larger than the 
common value of ai(Mjj) and a2(Mjj) (denoted 0:^2 (M*/)), then slightly larger 
values of my ,my will be allowed by the requirement of perturbativity for Yukawa 
couplings up to Mjj. When two- loop Yukawa couplings are included in the gauge 
evolution equations, an iterative procedure must be employed for finding a fully 
consistent solution. For a given mf ,rrib> choice and a given possible starting value 
of as{mz), full two- loop evolution of all three couplings may be performed (which 
requires evolving also the Yukawas at the same time). The value of Mjj at which a\ 
and «2 unify can then be determined. The common value, ai^iMjj), can then be 
compared to the evolved 0:3 (Mjj) value. This process is iterated until the evolved 
as(Mjj) value is as close to cxipiMjj) as possible without the Yukawas becoming 
too large (Aj < 3.3 is required). The precision of unification is then specified by the 
requirement that a^{Mjj) = ai^.{Mjj) to within a definite percentage deviation: 
for example unification to within 2.5% means that a^{Mjj) is allowed to be no more 
than 2.5% larger than oci^.{Mjj). Given experimental errors in sin 2 0w^ m z) an d 
variations in its best fit value as fourth-generation fermion masses and sparticle 
masses for all the generations are varied, we regard the region allowed by unification 
to within 2.5% as fully acceptable. However, to go much beyond the 2.5% allowed 
region of Fig. 2 would almost certainly require accepting the fact that one of the 
Yukawa couplings becomes non-perturbative at a scale below Mjj or that as(Mu) 
truly exceeds ai^iMjj), e.g. due to non-renormalizable operators, string threshold 
effects and/or evolution between Mjj and M atr i ng . 



I Not shown is ajj', the shift in au due to inclusion of two- loop Yukawa contributions is 
comparable to the a 3 (mz) shift, i.e. ^ 4%. 
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Of course, as at one loop, the value of 0:3 (mz) consistent with unification is 
sensitive to the effective scale Msusy implied by the sparticle masses. As the 
value for Msusy is raised, a^(mz) is reduced; for N g — 3 and Msusy — ITeV, 
the shift can be as large as ~ —10%. Similar results hold for N g = 4 as well. Given 
the present measurement of a^(mz) = 0.12 ± 0.01, it would seem that both the 
N g = 3 and N g — 4 scenarios predict an as (mz) which is somewhat high unless 
Msusy is significantly above mz- However, a large value of the effective M$usy 
is not easy to achieve, 1161 when the slepton, neutralino and chargino mass scales are 
significantly lower than the squark and gluino mass scales. Nonetheless, it is fair 
to say that 0:3 (m^) per se does not discriminate between the N g = 3 and N g — 4 
cases, since the increase in a^(mz) for N g — 4 vs. 3 from pure gauge effects can be 
compensated in the N g — 4 case by the large Yukawa effects in the gauge running, 
as just discussed. In addition, we shall find that N g = 4 scenarios characteristically 
force all the sleptons and squarks to higher masses (including those of the first three 
families), implying that the effective Msusy could be larger for N g = 4 than for 
N g = 3. Given this and the additional uncertainties associated with weak and 
GUT-scale threshold effects, we believe that the N g — 3 and N g — 4 scenarios are 
equally admissible. Attempts to address the moderately high prediction for az(mz) 
(that emerges in both cases) have recently been considered in Refs. [16,18]. 

In general, it is important to note that two-loop contributions to the running 
of a s (Q) will be much more significant in the case of N g = 4 than for N g = 3. This 
is simply due to the fact that the one-loop beta function for a s is proportional to 
(9 — 2N g ) (once Q > Msusy) and therefore is rather small for N g = 4. This means 
that two-loop contributions can represent a much larger percentage of the total beta 
function than in the N g — 3 case. The importance of two-loop contributions will 
be especially apparent when considering the running of the squark masses, where 
some terms involve the running value of a^. The cumulative effect of two-loop 
contributions can be large in such cases when evolving all the way down from 
Q = Mu to Q = mz- The situation for a\ and «2 is quite the opposite. Indeed, 
the one-loop beta function for 0:2 increases by a factor of 3 in going from N g — 3 
to N g = 4, implying that a,i increases much more rapidly as the energy scale is 
varied from mz up to Mu- 

What about Yukawa coupling constant unification? In Fig. 4, contours of 
constant \y(Mu), Xf(Mu) and A T >(Mu) within the allowed parameter space region 
are shown. From these contours we discover a number of important facts. First, 
Ay(Mu) and A T >(Mu) are generally quite different; Yukawa unification does not 
generally occur, although we see that in the my, my ~ 100 GeV corner of the 
tan/3 = 1.5 plot we do have Xf ~ \y ~ X r '- Second, we see that most of the 
non-perturbative borders are defined by one of the fourth-family Yukawa couplings 
becoming non-perturbative. At tan j3 = 1.5 the right-hand border results from non- 
perturbative behavior for \y (at the upper boundary, \ t becomes non-perturbative) 
while at tan/5 = 2.2 the right-hand border results when X T i becomes large and 
the upper boundary arises from non-perturbative behavior of Xf. However, it 
is also clear that, in general, not all of the fourth-family Yukawas (or At) are 
simultaneously large. This means that one is unlikely to be particularly close to a 
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Yukawa fixed-point solution' 211 such that the low-energy values of the A's are rather 
independent of their Mjj-scale values. The nearest approach to a fixed point occurs 
in the m? ~ my ~ 100 GeV corner of the tan/3 = 1.5 figure; the actual fixed point 

location is indicated by an X. In our earlier work) 11 with mt(mt) = 160 GeV and 
two-loop Yukawa contributions to gauge running not included in determining the 
parameter space, the large mf,my corner of the allowed region was nearer to this 
fixed point. 

3. EW Symmetry Breaking 

One of the many nice features of the MSSM extended by minimal SUGRA is 
that, for a very large region of soft-SUSY-breaking parameter space, radiative EW- 
symmetry breaking is automatically induced by renormalization group evolution. 
Thus, the hierarchy between Mjj (where scalar masses are universal and (therefore) 
EW symmetry is initially unbroken) and mz receives a natural explanation. In the 
three-family model, this EW-symmetry breaking is mainly a result of the quantum 
corrections arising from the large Xt Yukawa coupling which drives the i/2-Higgs 
field squared-mass to ever smaller values as the energy scale is decreased; EWSB 
occurs when the Higgs mass-squared term in the scalar Higgs potential is finally 
driven to a negative value. Of course, to obtain the precise value of mz as given 
by m 2 z = \{g 2 + g /2 )(vf +^2), there must be a relation among the supersymmetric 
model parameters. Mathematically, these relations result from requiring that the 
first derivatives of the scalar field potential with respect to v\ and v 2 vanish. Nor- 
mally, one of these conditions is used to determine the magnitude (but not the sign) 
of fj, in terms of the other initial parameters of the theory, including the soft-SUSY- 
breaking parameters. The other condition is used to determine the magnitude of 
7723, appearing in the -m\{HiH 2 + h.c.) mixing term in the scalar field potential. 
Thus, n and m\ at the weak scale become functions of mf, my, tan/5, mz, mo, 
m l j 2 and A. Since this radiative breaking mechanism is essential to the viability 
of the MSSM plus minimal SUGRA, we describe its dynamics for N g — 4 before 
turning to predictions for the SUSY sparticle spectra, and experimental constraints 
from direct particle searches and EW precision measurements. Not surprisingly, we 
find that in a four-family scenario the additional large Yukawa couplings also feed 
into the running of the Higgs field squared-masses, and modify the running signif- 
icantly. To indicate the way in which radiative breaking occurs in the four-family 
case, it is convenient to present the discussion at tree-level. A tree-level discussion 
is adequate for general understanding, and is a good approximation so long as the 

[22] 

tree-level minima obey all of the necessary stability and consistency constraints. 
However, we emphasize that our full numerical calculations are actually performed 
by minimizing the full scalar potential at one-loop. 

Figures 5a and 5b illustrate how radiative breaking occurs in the tree-level 
approximation in one specific four-family case. We take mt = 165 GeV, tan (3 = 1.5, 
mf = 100 GeV, my = 100 GeV and m T ' = m v i = 50 GeV. As outlined in the 
introduction, we specify our soft-SUSY-breaking parameter boundary conditions 
in terms of £0 = mo/my 2 , £a = A/my2 and m^- For this illustration we 
have adopted the string-motivated dilaton boundary condition scenario specified 
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by £o = l/"\/3, £a = ~~ 1- I 11 the plots, m 2 Hv m 2 H2 are the Higgs mass-squared 
parameters appearing in the soft-SUSY-breaking potential (at Mjj they are equal 
to ttiq), while m 2 2 — m % 12 + A* 2 are the parameters multiplying the iff and iff 
quadratic scalar field terms in the full scalar field potential; as defined earlier, 
m 2 , is the coefficient of the H1H2 mixing term in the scalar field potential. In 
Fig. 5, we see that the mjj. start out above zero at Q — Mjj, and then evolve 
below Mjj so that both eventually take on negative values (but with rn? H2 being the 

more negative). The \i parameter is determined at the weak scale by minimization 
of the RGE-improved tree-level Higgs potential, and is then evolved up to Mjj. 
The fact that (J?(mz) > for the chosen value of tan/3 and the given value of 
mz indicates that an allowed four-family EW symmetry breaking solution via the 
radiative breaking mechanism exists. 

However, some differences in comparison to N g = 3 are apparent. Even for 
the small tan/3 = 1.5 value, m\ x has significant evolution due to the large A&/ 
and X T i Yukawa couplings. In the standard three-family case, m\ x only evolves 
significantly when tan/3 is so large as to require a large value for A5. In Fig. 5b, 
one can see that although the scalar potential parameters m 2 , m\ and m\ each 
evolves separately and ends with a value that is > 0, in combination they serve 
to trigger the breaking of the EW symmetry, measured by the tree-level stability 
condition S = m\m\ — m\. It is of course possible for the mechanism to fail, i.e. 
if /x 2 < 0. This can happen in the case where m\^ is too large compared to m^. 

4. The N g = 4 Higgs Sector 

We begin* by reminding the reader that at tree-level the Higgs sector is deter- 
mined by just two parameters, tan /3 and m^o (the mass of the CP-odd scalar Higgs 
boson), and that if m^o is large (as we shall see it is in all the SUGRA models 
considered) then the lightest CP-even Higgs boson hP is very SM-like and has mass 
bounded from above by mz- However, it is well-known that one-loop radiative cor- 
rections to m/jo can significantly increase the upper limit of m^o (denoted m™, ax ) 
for the large mt value found in the CDF and DO experiments. Exactly how large 
the upper bound is depends upon other SUSY parameters, the most important sen- 
sitivity being to the stop squark mass, m-r. In the absence of a fourth generation, 
if the running masses are mr — 1 TeV and mt = 165 GeV, then m^o can be as large 
as 120 — 125 GeV at large tan/3. However, in many SUGRA models (for example 
the dilaton scenario) mo and, consequently, m^-are significantly smaller than 1 TeV 
for m~ values below 1 TeV. In such models m™ ax is typically ^ 100 GeV. 



* An excellent brief overview of Higgs phenomenology and discovery techniques is now avail- 
able in the Higgs subgroup summary appearing in Rcf . [23] . A longer version of this review 
will soon be available, Rcf. [24]. References for statements not explicitly referenced below 
can be found in these reports. 
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The magnitude of m A o determines the observability of the A and H° Higgs 
bosons, which are more or less degenerate when m A o is large. For example, e + e~ 
colliders can only probe up to m A o ~ m H o ~ ^/2 - 30 GeV in the Z* -> H°A° 
mode (the only viable mode when hP is SM-like). For N g = 3, SUGRA model 
predictions for m A o can range from 200 GeV on up, and an e + e~ collider with 

= 500 GeV would at least have a small chance of seeing H° + A production. 
The parameter \x is also of interest in that its origin in the SUGRA models is rather 

uncertain. For N g — 3, pL tends to take on rather moderate values ~ 500 GeV. 1161 
We now discuss the influence of a fourth generation. 

First, we note that a fourth family will give an additional set of t',t', v 1 ,P ', 
b', b' and t', t' loop contributions to the radiative corrections for m™f x . (As noted 
earlier, as a consequence of the many loop corrections for N g = 4 it is essential 
that the one-loop effective potential be used in the scalar potential minimization 
as well as in the determination of the Higgs masses.) A fourth family will also 
typically lead to a rather high m~ scale, as discussed in the following sections; this 
will in turn influence m^o and //, which will often take on relatively large values. 

Typical results for /i, m A ° and m^o are illustrated in Fig. 6, where we display 
contours of constant /i, m A ° and m^o in my, my parameter space for m t (m t ) = 
165 GeV, tan/3 = 1.5 and = 600 GeV in the dilaton scenario: £o = I/a/3, 

£ A = —1. We see that /x is generally above 1 TeV, a substantial increase over the 
corresponding N g = 3 dilaton scenario result; m A o is typically also rather large, 

ranging from ~ 600 GeV to above 1 TeV — a yHs = 500 GeV e + e~ collider would 
not allow detection of the Z* — > H°A° pair production process. The large values of 
m A o ~ mjjo imply that one must also be cautious to account for supersymmetric 

decays of the A and H°. These decays can deplete the more easily observed bb 
decay channel, especially given that tan/3 cannot be large for N g = 4, implying 

that the H°, A — > bb coupling cannot be greatly enhanced relative to the SM-like 
result. We will not attempt a detailed study of the H° and A decays here. 

The above results are not greatly altered by changing the Mjj boundary con- 
ditions, keeping m^ 2 fixed at 600 GeV. In the high-mo (Co = 1) scenario, both 
m^o and /1 become somewhat (~ 100 GeV) larger over most of the my, my param- 
eter space. The only exception is the corner where my is big and mf is small; in 
this corner m^o is slightly smaller (564 GeV compared to 610 GeV) in the high- 
mo scenario than in the dilaton scenario. In the case of the no-scale boundary 
conditions, £0 — £,A = 0, both \x and m^o move to lower values (shifts are of 
order 50 — 100 GeV). For example, in the small mt>, large my corner one finds 
m A o ~ 450 GeV. 

Fortunately as regards the prospects for h° detection, the extra loops from 
the fourth family do not yield overwhelmingly large radiative corrections to m™f x . 
They increase m^o ax by 10 to 25 GeV relative to corresponding N g = 3 predic- 
tions. As illustrated by the contours in Fig. 6, at worst m^f x ~ 130 GeV at the 
perturbative boundary in the my,m t > parameter space. For all the three Mjj sce- 
narios, m/jo remains very much in the 115 — 122 GeV range for most of allowed 
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my, my parameter space, for this relatively large m^ 2 = 600 GeV value. (Lower 
m/jo values are predicted at lower m^ 2 -) The phenomenology of the h° depends 
upon whether or not it is SM-like. For the bulk of parameter space, and certainly 
for the preferred scenarios and portions of parameter space, m^o is large and the 

h° will be very SM-like. We discuss its phenomenology assuming that this is the 
case. 

Because m^o is predicted to be ^ 100 GeV for many N g = 4 scenarios, LEP 
II would be less likely to find the hP if there is a fourth family. However, a y/s = 
500 GeV e + e~ collider would have no difficulty in doing so in the Z* — > Zh° 
production mode. At the TeV* upgrade of the Tevatron, detection of a SM-like 
h° is probably only possible in the Wh° — > ivbb mode, and then only if m^o < 
95 GeV. (Although it is not impossible that the Wh° — > £vt + t~ mode could be 
used for 110 ^ m^o ^ 120GeV. [2jl ) For the many four-family scenarios that lead 

to m^o above 100 GeV, searches for the h° at the Tevatron would be, at best, 
problematical. 

At the LHC, a SM-like Higgs boson (for which we use the generic notation 
h below — the hP might or might not be perfectly SM-like) in the mass region 
being discussed would typically be found through production via gg — > h and 
decay to either 77 or ZZ* (with ZZ* — * 4£). The gg — > h production rate, 
proportional to F(h — > gg) would be greatly enhanced by the additional t' and 
b' loop contributions to the one-loop gg — > h coupling. This is illustrated in 
Fig. 7 where we plot for a SM-like h the ratio of N g = 4 to N g = 3 values for: 
T(h — > gg), T(h — > 77), and F(h — > gg) x BR{h — > 77), taking tan/5 = 1.5, 
m t = 165 GeV, m t > = my = 100 GeV, m T > = m v i = 50 GeV, and assuming 
that superpartners are sufficiently heavy that their contributions to these one-loop 
quantities are small. (As discussed in Sec. 7a, experimental limits tend to prefer 
this type of scenario.) For F(h — > gg), and hence the gg — > h production rate, we 
see an enhancement by a factor of 10. Thus, the 4£ channel, which for N g = 3 is 
only viable for m^ f& 130 GeV, would yield a detectable signal down to somewhat 
lower masses, perhaps as low as m^ = 120 GeV. (BR(h — > ZZ*) falls very rapidly 
with decreasing m^ so very few events would result for m^ values much below this.) 

The h — > 77 decay also arises at one-loop. For N g = 3, the main contribution is 
from the VF-loop diagram. Fermion loops (for massive fermions) cancel against the 
VF-loop contribution and decrease the h — > 77 width. For N g = 4 this cancellation 

can be quite substantial, as illustrated in Fig. 7 [ . 26] Combining the resulting reduc- 
tion in BR(h — * 77) with the enhanced gg — > h production rate, the resulting 77 
channel event rate tends to be substantially suppressed relative to the N g = 3 rate 
in the 100 — 130 GeV mass range of interest. In fact, the preferred m?, my, m T ' 
mass choices delineated above are about the worst that can be made in this regard. 
Thus, detection of the h° in the 77 channel at the LHC becomes problematical, 
even when m^o is large and the h° is SM-like. 

Not investigated to date is whether the enhanced production rate from gg — > h 
might make detection of a SM-like h° in the inclusive bb channel possible (assuming 
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high 6-tagging efficiency and purity). Finally, we note that the extra family has 
essentially no impact on the pp — > tih° — > ibbbbX LHC detection mode, which 
would continue to be viable for m h o < 120 GeV, when the hP is SM-like. [231 

Regarding the heavier A and H°, for our typical scenario they are so massive 
that perhaps the only accelerator with adequate energy for their production and 
possible detection will be the LHC. At the LHC, for N g = 3 the detection of a 

massive A or H° is possible only if tan (3 is so large {e.g. <^ 10), that the gg — > bbH° 
and gg — > bbA° production rates significantly exceed the inclusive gg — > _£/~°,A° 
rates! 23 ' 241 However, for N g — A the 6' and t' loop contributions to gg — > A and 
— > if will greatly increase these inclusive production rates, regardless of the 
tan (3 value. The possibility of observing the H° and A at the LHC in the inclusive 
bb final state assuming such highly enhanced rates should be carefully examined. 

[27] 

As a final aside, we note that the standard relation, 



2 rnl + fi 2 
sim p 



valid for tan/5 & 20 if N g = 3, is not necessarily maintained for N g = 4 since there 
will be corrections involving the \y,\ T ',\b Yukawas (where A& can essentially be 
neglected). In Fig. 6, for mt'^my large, e.g. m t > = my = 100 GeV, the relation is 
satisfied to within 3%; however for my ~ 110 GeV, my = 55 GeV, the relation is 
seriously violated: {m? A0 sin 2 /3)/(m~ + /i 2 ) ~ 0.25. This latter situation is realized 
in the region where my > my, a region that is experimentally disfavored unless an 
unnatural quark mixing pattern exists (see Ref. [1] and Sec. 7 for a more detailed 
discussion of this point). 



5. A Fourth-Generation Sparticle Constraint on the SUSY Scale 

We turn now to an important additional constraint on soft-SUSY-breaking pa- 
rameters that can arise from consistency of the sparticle mass spectrum with LEP 
limits and a neutral LSP (lightest supersymmetric particle). (A charged LSP is 
excluded experimentally.) A glance at Eqs. (A. 11) through (A. 19) in the Appendix 
shows that dm^/dt (where t = ^ ln[Q( GeV)]) for squarks and sleptons receives 
positive contributions from Yukawa terms and negative contributions from gauge 
terms. Thus, starting from a universal m\ and evolving downwards in t to mz, 
the lightest squark/slepton will be the one with the largest Yukawa contributions 
relative to gauge contributions. This turns out to always be the t' r or the v' R . After 
including r' R — t' l mixing, the lightest r' eigenstate is denoted t[, and similarly u[ 
is the lightest V eigenstate. 

The u[ can be even lighter than the r[ in scenarios with small £o, such as 
the dilaton and no-scale models. However, we see no general phenomenological 
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reason for not allowing the u[ to be the LSP. Indeed, it is even quite likely that 
the Xi decays invisibly when m~o > m~/ via ~^ v r^'\-, assuming at least a small 
non-zero value for the required 3 — 4 generation mixing angle. Thus, we will only 
impose a phenomenological limit on the charged t[. 

Eq. (A. 15) shows that the Yukawa contributions to dm 2 , Idt (in the Appendix 

T R 

we use the notation t' r = E') are controlled by V T > = m 2 Hi + m~, + m~, + A 2 ,. For 
moderate initial values, m~^ = m~^ = mo, the more negative m 2 Hi is (see Fig. 5) 
the larger will be m~, ; m 2 H in turn becomes more negative for larger my. Thus, 

R 

for some choices of soft-SUSY-breaking parameters it is possible that m~/ will fall 
below the 45 GeV LEP limit and/or below the mass of the lightest supersymmetric 
particle (the lightest neutralino, Xi, m the models we consider) for low values of 
my (but not for higher values). 

To illustrate this, we return to Fig. 6, where we have chosen tan/3 = 1.5, 
mt(mt) = 165 GeV and dilaton boundary conditions for the soft-SUSY-breaking 
parameters with = 600 GeV. The diamonds indicate the portion of parameter 
space at low my that is ruled out because one predicts m~/ < 45 GeV or m~/ < m~o. 
(For this particular choice of m^j m x\ ~ 108 GeV and it is the m~/ > m~o 
requirement that fixes the diamond region.) By raising m^rt, this problem region 
is moved to lower my values since the starting value of m\ oc m^ 2 for m~, is 
increased more rapidly than the off-diagonal mixing term A oc m^ 2 - Conversely, 
by lowering m^ 2 we eventually reach a value for which no portion of parameter 
space remains allowed. 

The type of boundary condition applied is also important. For instance, for the 
same = 600 GeV but £o — 1, £,A — ~ 1 (the high-mo scenario) no points with 
my,mf > 50 GeV are eliminated by virtue of the r{ constraint, whereas for the 
no-scale choice of £o — £,A — 0, the portion of parameter space removed expands 
to include slightly higher values of my than in the dilaton scenario. 

Let us further expand upon this point. From the above discussion we see that 
in the generic one-loop formula, mf = m\ + (7jm^, 2 + Dim 2 z cos 2/5, the evolution 

parameter C{ for % — t' r is strongly affected by indirect effects from Ay, due to the 
absence of an a s contribution to C{. In general, larger my tends to raise C~> . As wc 
have discussed above, this can be traced to an increasingly negative contribution 
in the RGE for m~, from the increasingly negative value of m 2 H as my increases. 

Fig. 8 demonstrates this sensitivity of the m^-scale value of m~, to my in the case 
of m t > = 100 GeV, tan/9 = 1.5, m 1 / 2 = 600 GeV with dilaton boundary conditions: 
as my is lowered from my = 110 GeV to my = 90 GeV, the m^-scale value of 
m~, is lowered by ~ 40%, with the lowest r' physical eigenstate mass decreasing 

T R 

dramatically, from m~, = 311 GeV to 189 GeV. In contrast, the evolution and 
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m^-scale values of and m~, are little affected. Note the subtlety of the log(Q) 
behavior of m~, . In the region of small (and decreasing) \og(Q) m~, rises due 

T R T R 

to the increasingly more negative value for m 2 Hi as \og(Q) decreases (see Fig. 5). 
Since this rise is less for smaller my, for small enough my one can even obtain 
m~, < 0, and EM will no longer be unbroken. However, as my is decreased one 

T R 

first arrives at a point where either m~/ < m~o and the LSP is no longer neutral, 

or m~/ < 45 GeV, violating LEP limits. We re-emphasize that these requirements 

result in the strongest additional RGE-related phenomenological constraint deriving 
from the SUGRA extension of the four-family model. 

To gain some additional insight regarding the stringency of this bound as a 
function of M(y-scale boundary conditions, we present the right-hand boundaries, 
analogous to that of the diamond region shown in Fig. 6, for a selection of different 
possibilities in Fig. 9. The right-hand window shows results for the dilaton case of 
(f , £ A ) = -1) for m 1/2 = 900, 600, 300, 250 GeV. Very little of the Yukawa- 

allowed parameter space survives for the lowest choice. In the left-hand window 
we fix m\j2 = 400 GeV and vary (£o,£a)- For the no-scale choice of (0,0) very 
little of parameter space yields an acceptable r{. 

Of course, the m~^ > max |45 GeV, m~o| constraint is significant even if we 

adopt a large (but allowed) value of my and do not employ a specific scenario 
for £o an d £a- This is illustrated in Fig. 10. There, we take (m(m) masses) 
m t = 165 GeV, mf = my = 100 GeV, m r > = m v > = 50 GeV and tan/5 = 1.5 (for 
these mass choices, solutions are only allowed for tan/5 values very near 1.5) and 
plot the smallest possible value for £o = mo/ m i/2 that is allowed as a function 
of mi/2 after scanning over £4 in the range [—3, +3]. A very important generic 
feature emerges from the steep rise of the minimum £0 value as m^i^ decreases: 
for given my , my , m T > , m v i there is a definite lower bound on m\j2 (for reasonable 

values of £0) arising from the m~/ > max |45 GeV, m~o| constraint. For the 

t' and b' masses considered, this bound is m\i2 <S 140 GeV. If we recall that 
m~ ~ (as(mz) /ajj)m 1 / 2 , and that a^(mz) / 'oljj ~ 1.4, we see that this m~ bound 
translates into a significant lower bound of m~ ^ 200 GeV deriving purely from 

limits on the fourth-generation t[ mass. This type of constraint does not arise in 
the N g = 3 MSSM. 

The ultimate lower bound on m\j2 (along with the corresponding lower bound 
on m~) is actually quite independent of my. For instance, if my is lowered to 
50 GeV, keeping mt and mp fixed, at £0 — 5 the lowest allowed value of 

decreases by only about 5 GeV compared to the my = 100 GeV value of ~ 140 GeV 
illustrated in Fig. 10. This is because, at the ultimate lower bound, mo is very 
large and the m 2 Hi term in U T > (see above) is swamped by the m~, and m~, terms 

which start off of order m\ and remain large. 
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We re-emphasize the fact that in a specific SUGRA scenario with a fixed value 
of £o (and of £4) the lower bound on m^ 2 (and m~) can be much larger than the 
high-^o ultimate lower bound if £0 is small, as is illustrated in Fig. 10. (Of course, 
whatever the lower bound on m~{m~), m~{pole) will be roughly 5-6% higher.) We 
also note that at fixed £0 there is some dependence of the lower bound on roy 2 
upon the value of £4. This dependence is extremely weak at high £0, where £4 
values in the entire [—3, +3] range generally give an allowed solution for the lowest 
acceptable value. But for lower £0 values, the lower bound on my2 is generally 
achieved only for £4 values near 0. (In fact, the appropriate £4 range more or less 
scales with the magnitude of £)•) 

Perhaps one other plot is useful in fully understanding the t[ constraint. In 
Fig. 11 we consider various quantities as a function of £4 (ranging from —3 to +3) 
for fixed values of £0 i n the range from to 5. In all cases we take my (my) = 
Tnt'(mt') = 100 GeV (corresponding to pole masses of approximately 105 GeV). 
For each £},£4 choice, we determine the minimum value of my2 (the scan being 
confined to the region m^ 2 < lOTeV) that is consistent with the r{ constraint. 
We plot the (minimum) value of m~{m~) corresponding to this minimum value of 
my 2- (F° r later phenomenological use, we also plot the corresponding m^o and 
my {pole) — m/jo — values.) For curves of limited extent in £4, the termination 
point (s) define the range beyond which consistent solutions are not found with 
mi/2 < lOTeV. The inconsistencies that arise at large £4 are of two types: (i) 

that the EWSB solution requires m 2 A0 < 0; and/or (ii) that a 4th-family sparticle 
with color or charge must have m 2 < 0, thereby breaking the color and/or U(l)- 
electromagnetic symmetries. 

This plot defines precisely the minimum m~ value that one can have for a 
given £oj£a boundary condition choice. Once again we see an absolute lower limit 
of order m~{m~) ^ 200 GeV for our choice of my,mt>] the value of tan j3 for which 
the minimum m~ value is reached is always in the range 1.5 to 1.6. Note that for 
low £), the allowed solution range for £4 is limited, and that the minimum m~ 
achievable increases substantially when £4 is not near 0. 

6. The Sparticle Spectrum 

We shall find in Sec. 7 that the first evidence for a fourth family is very likely 
to be discovery of the b', t' ', r' and/or v' at LEP-II and/or the Tevatron. If 
one or more of these fourth-family members are found, the immediate question 
will be how this impacts the supersymmetric particle spectrum, especially in the 
standard renormalization group equation (RGE) context. It is this latter issue that 
we address in this section. First, we highlight the main features of the sparticle 
spectrum in a four-family scenario, and compare the results to those obtained in 
the three-family case. As previewed in Sec. 4, we find that it is more than likely 
that some of the fourth generation squarks and/or sleptons will be lighter than 
their counterparts in the first three generations, with the t[ most probably being 
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the lightest. However, testing consistency of their masses with unification and the 
RGE's is likely to be challenging given the possibly large A-term induced mixing, 
and the probability that they will have strange decay patterns and be tricky to 
observe. Thus, we will focus primarily on gaugino masses and on the masses 
of the sleptons and squarks of the first two generations, focusing on when and 
how correlations among these masses will be indirectly sensitive to the presence 
of a fourth generation through the RGE's. Certain relationships between masses 
are rather insensitive to whether N g = 4 or N g = 3, and thus provide a test 
of the general RGE context and universality of boundary conditions, while other 
mass correlations are very different depending upon the value of N g . These latter 
relationships with strong N g dependence would provide indirect evidence for the 
presence of a fourth generation, even if no particle or sparticle belonging to the 
fourth generation is directly observed. 

We first present some sample mass spectra, then discuss mass sum rules and 
relations, and finally focus on a specific correlation between the first-family slepton 
masses and the LSP mass that could reveal the presence of a fourth generation. 



6a. Sample Mass Spectra 

For our illustrations we shall adopt \i > 0, mt> = my = 100 GeV, and m T > = 
m v i = 50 GeV. We scan over the allowed tan /3 values at any given choice for 
m~(m~) and plot mass spectra in units of m~(m~) . Three Mjj boundary conditions 
will be considered: 

1. the dilaton scenario, with £4 = A/m 1 / 2 = — 1, £0 = m o/ m i/2 = Va/3; 

2. the no-scale scenario, with £4 = £0 = 0; and 

3. the high-mo scenario, with £4 = —1, £0 = 1. 

We begin by focusing on the spectrum for the r[ for these three models. 

In Fig. 12, we plot the m~^ spectra obtained by scanning over allowed tan/5 
values at various m~ values. The lower limit on m~> is fixed by m~> > 45 GeV in 
the case of the high-mo scenario, and by m~^ > m~o in the no-scale and dilaton 
models. The dramatic decrease of m~, as m~ decreases, discussed in Sec. 5, is 
evident. Of course, at high m~ values, the value of m~^ is given roughly by m~, ~ 

m\ + C~,m\i 2 oc (£q + C~r)m~, which is larger for larger values of £o- 

Let us now turn to the typical mass spectra for all the other supersymmetric 
particles. We begin with results for dilaton boundary conditions: Figure 13 shows 
our results. These can be compared to the N g — 3 results from Ref. [5] for this 
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same boundary condition choice. 

The first big difference between the N g = 4 and N g — 3 cases is the much larger 
lower bound on m~ in the former case. As discussed previously, this reflects the 

m~/ > max |45 GeV, m~o j requirement that is violated for low m~ values. A less 

dramatic difference is that the allowed range of tan (3 is smaller for N g = 4, and 
thus the value of m~ almost completely fixes the chargino and neutralino masses, 
whereas for N g — 3 there is some scatter. The limitation of tan/3 to low values 
for N g = 4 also means that (unlike for N g = 3) the t\ remains almost degenerate 

with the £rs of the first two families; thus, a separate plot for the t\ is not given 
in our figures. At a given m~, the most significant N g = 4 vs. N g = 3 difference is 
the much larger masses for the squarks and sleptons of the first three generations. 
Finally, there is the simple fact that fourth-generation squarks and sleptons are 
present for N g = 4. We have already noted that t[ and u[ tend to be the lightest 
of the sleptons. From Fig. 13 we see that the lightest squark is very likely to be 
the b[. This is a rather general result. 

As already hinted and more directly demonstrated in the next subsection, 
the large squark and slepton masses for the first three generation members can 
be directly traced to the much larger value of a(Mjj) for N g = 4. The basic 

idea is that mf = Wq + QmL + -D;m| cos 2/3, where i indicates the squark or 

slepton in question and m ly / 2 = m g a {Mu) / ' 0:3(777,-) . Since the C\ are not terribly 
different in the three- and four-generation cases, the much larger a{Mu) in the 
N g = 4 case greatly increases the squark or slepton mass at given m~. Combining 
this effect with the higher lower bound for m~ in the N g = 4 vs. N g — 3 case 
results in a large increase in the lower bounds on squark and slepton masses in 
going to N g = 4. For example, mj^m~,mj R ^ 250, 240, 200 GeV compared to 

mj^,m~, mj R <; 110, 52, 45 GeV in the N g = 4 and N g = 3 cases, respectively. The 
high lower bound values in the N g = 4 case place these sleptons from the first 
two families beyond the reach of LEP-II and almost beyond the reach of a ^J~s = 
500 GeV NLC. We shall see that the large squark and slepton masses relative to 
gaugino masses might well provide the most compelling indirect indication for the 
presence of four generations that one can obtain using only particles and sparticles 
belonging to the first three generations. 

In Fig. 14 we give the corresponding results assuming no-scale boundary con- 
ditions: £a = £o = 0. The first noteworthy point is the much higher lower bound 
on m~ that arises in the no-scale compared to dilaton boundary condi- 

tions, when N g = 4. Due to the zero value of mo at Mjj, the slepton masses 



* In making comparisons, it is necessary to note that for N g — 3 the value of — 0.12 was 
employed in Ref. [5], as compared to ~ 0.128 or so for our full two-loop four-generation 
treatment here. Since M\ : M2 : M3 ~ ol\ : a-i : olz (all at mz), rn~ , m~ , m~+ masses at a 

Xi X2 Xi 

given m~ in Ref. [5] for N g — 3 are approximately 0.94 <~ 0.12/0.128 times those appearing 
in Fig. 13. 
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are smaller than for dilaton boundary conditions, and, in particular, m~/ is more 
easily driven to too low a value. In addition, for the lower m~ values allowed by 
the m~/ > m~o constraint, the u[ can be the LSP; it is often very substantially 

lighter than the x^l- However, the and v masses for the first two families are 
still much bigger than found in the comparison results for N g = 3, see Ref. [5]. 

The In is significantly lighter than the £l in the no-scale scenario (whatever N g ) 
because the associated soft-mass-squared evolution is fed only by the U{1) gaug- 
ino mass terms, and moves to positive values (starting at Mjj from for £o = 0) 

much more slowly than the soft-mass- squared that contributes to the £l and u 
masses, which is fed by SU{2) as well as U{1) gaugino mass terms (see Sec. 6b). 
Squark masses remain very similar to the dilaton scenario results due to the fact 
that the C~ terms dominate the mo terms for strongly-interacting sparticles. As 
in the dilaton case, the large squark and slepton masses relative to m~ provide a 
signal for N g = 4. Again there is a large difference between N g = 4 and N g = 3 
in the lower bounds so crucial for LEP and NLC phenomenology. For example, 
mj L ,mi;, m j R ~ 265, 260, 150 GeV compared to mj^m-, rnj R ~ HO, 60, 70 GeV for 

N g = 4 and N g = 3, respectively. Thus, for N g = 4 only the £r would be within 
the kinematical reach of a ^/s = 500 GeV NLC. 

In Fig. 15 we give results assuming the high-mo boundary condition case: 
£o = — C.A = 1- Here, the r{ constraint is more easily satisfied and m~ can take 
on lower values. Even at the lowest allowed m~ values, the u[ is heavier than 
the t[ (see Fig. 12), which in turn is required to be heavier than the LSP; in 
fact, even at low m~ the u[ (but not the r{) is heavier than the xf . Squarks 
and sleptons receive a significant fraction of their mass from the large mo value, 
and are thus more similar in mass. However, the C~ contributions are still very 

important, and squark and slepton masses for members of the first two families 
continue to be much larger for N g — 4 than for N g = 3. To repeat our previous 
sample comparisons, mj , mj;, mj ^ 245, 240, 220 GeV for N g = 4 compared to 
mj L ,m~, mj R ^ 120, 89, 107 GeV for N g = 3. 

6b. Sum Rule and Sparticle Spectrum Tests for N g = 3 vs. N g = 4 

In the models considered here, having universal soft-SUSY-breaking scalar 
mass mo, there are many sum rules relating the gaugino, squark and slepton masses. 
Here, we survey the impact of a fourth family on the masses and mass sum rules 
for the gauginos and for the squarks and sleptons belonging to the first two fami- 

[28] 

lies; we also point out several differences between the sum rules that relate only 
fourth-family masses and those that involve only the third-family masses. Gener- 
ally speaking, because of the large mixings that can be present in both third and 
fourth generation squark and slepton mass matrices, the most precise tests of the 
consistency of RGE evolution and grand unification with sparticle mass spectra 
may be those employing first and second generation members, for which the mass 
matrices are very nearly diagonal. 
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To proceed, we must first discuss gaugino masses in more detail. We employ the 
notation k — 1, 2, 3 for the U(l), 577(2), and 577(3) groups; Qk{t) are the associated 
running coupling constants, and M k (t) are the associated running gaugino masses, 
where t — ^- log[Q( GeV)] - tjj denotes t at Q — Mjj, tz denotes t at Q — mz, 
and so forth. If the gaugino masses do indeed take a universal value m^ 2 at Mjj 
then at one-loop we have the standard relations: 



M k (t) = [a k (t)/a(tu)]my 2 ; 



(1) 



see Eqs. (A. 28) through (A. 31) in Appendix A. For an approximate idea of the 
numerics for N g = 4 as compared to N g = 3, we take a~^(mz) — 128, sin 2 9\y = 
.2316 and a 3 (m z ) = 0.132 at t = t z and a(t v ) = 0.0917 (0.043) for N g = 4 
(N g = 3). This gives: 



(Mi(t z )\ 

M 2 (t z ) 
\M 3 (tz)J 



JV =4 



/0.185\ 

0.37 
V 1-44 / 



mi/2 ; 



fM 1 (t z )\ /0.39\ 

iV =3 



M 2 (t z ) 
\M 3 (t z )J 



0.78 
V 3 / 



mi/2 • (2) 



These results are slightly modified when two-loop Yukawa contributions to the 
running of the gauge coupling couplings are incorporated. One finds coefficients for 
N g — 4 of roughly 0.194, 0.386 and 1.47 for k — 1, 2 and 3, respectively. The main 
effect derives from the simple fact that a(tjj) decreases so that the eti(tz)/cx(tu) 
ratios increase. For example, for = my = 100 GeV, a(tjj) is shifted down by 
about 5% to about 0.0874. This results in an increase of the N g — 4, k — 1, 2 
coefficients in Eq. (2) by about 5%; for k — 3, a 3 (tz) also decreases (to about 
0.1288) and the net Eq. (2) coefficient increase is only about 2%. 

We also recall that simple asymptotic results for the low-energy masses of the 
gauginos Xi an d xt arise if mw ± M 2 \ (but not if mw ~ \\f^\ ± M 2 \). 

These results are summarized in the chargino mass formulas and neutralino mass 
formulas given in Eqs. (36) and (37) (for charginos) and (40), (41), (42) and (43) 
(for neutralinos) of Ref. [29]. (See also Eqs. (5.3) and (6.8) of Ref. [28].) An 



example of such asymptotic results are the large a relations m~o ~ M\ and m~+ ~ 

A.1 Al 

m^o ~ M 2 , where M\ i2 are evaluated at energy scales of order mz- As seen from 
Figs. 13, 14 and 15, the above approximate mass formulae work reasonably well, 
but not perfectly, at large m~. For instance, at high m~ one finds m~o/m~ ~ 0.135 

compared to M1/M3 ~ 0.131 (with m? = my = 100 GeV two-loop effects included 
in the gauge running). 

One of the first priorities of a next linear e + e~ collider (NLC) and the LHC 
will be to test the M^j-scale universality assumption for the Mj's. At the NLC it is 



estimated that m 



m~o and m~o can be measured to within a few GeV 

Al A2 



[30] 



This 
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will allow a pretty good determination of Mi(tjj) / M 2 (tjj) ■ The optimal situation 
arises if // is large, as is essentially always the case in three-generation models, and 
very often the case in the four-generation models. Then, as summarized above (and 

are of order mz) m~o ~ M\{tz) and 



making the approximation that m~o and m~+ 

Al Al 

m~+ ~ m^o ~ M 2 {tz) and we have the more general result: 



Mijtu) _ a 2 (tz)M 1 (t z ) _ 
M 2 (tu) a 1 (t z )M 2 (t z ) ~ ai(t z )m~4 



a 2 (t z ) m x\ 



(3) 



(For greater accuracy, one will wish to include the two-loop corrections to this 
relation when actually performing this test.) Our ability to test the universality 
for % = 3 will probably be much more limited. At the LHC, m~ is unlikely to be 

measured to better than 50 to 100 GeV. Further, due to the small beta function 
at one-loop for i = 3 in the case of four generations, two-loop corrections to the 
running of 0:3 can be significant, and these depend at least somewhat upon other 
model parameters. Nonetheless, it should be possible to extract a ±(20 — 30)% 
value for M 2 (Mjj) / 'M^(Mu) . The determination of the Miitjj) is crucial, not only 
as a test of universality, but also in making predictions for the squark and slepton 
masses and testing mo universality. 

But before turning to squarks and sleptons we must make a few more comments 
on the gaugino sector. First, the universality test will not have any significant 
sensitivity to N g = 4 vs. N g = 3. Second, although the absolute mass scales are 
sensitive to the limits on m\j 2 discussed earlier that keep m\j 2 from being as small 
for N g = 4 as it can be for N g = 3, it is difficult to use the absolute mass scale 
for a reliable probe of the value of N g since N g = 3 models with N g = 4 type re- 
values are certainly entirely viable. Sum rules such as 



m~ + m~ t = \M 2[ x 



m\ sin2/3j , 



(4) 



following from the determinant of the mass matrix, are sensitive to the N g = 4 
restriction 1 ^ tan/3 ^ 3 if M 2 /j, ~ rriyy But such small values are not the norm, 
and to use this particular sum rule we must detect both xf and x 2 ■ Thus, we must 
include squark and slepton masses in our considerations in order to gain sensitivity 

tO Ng. 

We shall see below that the squark and slepton masses at one-loop are deter- 
mined by the mo soft scalar mass (or masses should the tuq values not be universal), 
by (5 and by the functions: 



C k {t)m\ /2 = 



il\ 



VI/ 



tu 



X 



7T 



dtgl{t)M 2 k {t) 



(5) 



If one-loop evolution is used, then the integrals can be analytically carried out to 
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give: 



C k {t)m\ /2 = 2 



fi\ 

3 
4 

\v 



M 2 k {t) 



1 - 



a 2 (t 



v , 



«*(*) 



/ 9-2iV ff 



with 6^ 



5-27V 

3 



9 

V-|-2^. 



(6) 



where the are the one-loop beta function coefficients for the gauge couplings: 
da>k/dt = —b k a\. If we input the relations of Eq. (1), then (at one- loop) we find: 



Cfc(f) 



a 2 (tu) 



K 1 



(7) 



However, for four generations, one-loop evolution for 0:3 is generally not an ade- 
quate approximation when calculating C3; the evolution of 0:3 is sensitive to two- 
loop terms, due to the small value of 63, and a 3 appears to the third power in the 
expression for C3 in Eq. (5). The values of for N g = 4, computed numerically 
with full two-loop evolution for the couplings, are compared to those for N g = 3 
in Table 1 at scales mz and 1 TeV. For the N g = 4 computations we adopted 
m? = my — 90 GeV, m T i = m v i = 50 GeV, and tan/5 = 1.5. Results are fairly 
insensitive to these choices. In the case of N g = 4, the results for Ck would have 
been roughly 20% bigger had we neglected two-loop Yukawa terms in 0:3, which 
keep 0:3 somewhat smaller (as described earlier) than otherwise. 



Table 1: We tabulate the Ck values for N g = 4 and N g — 3 at the mz and ITeV energy scales. 
For N g = 4, we employ m t i = my = 90 GeV, m T t = m v , = 50 GeV, and tan/3 = 1.5. 



k 


C k (tz);N g = A 


Cfc(tiTcv); N g = 4 


C k (tz);N g = 3 


Cfc(^lTev); N g — 3 


1 


0.131 


0.130 


0.151 


0.148 


2 


0.376 


0.371 


0.484 


0.459 


3 


4.59 


4.06 


7.30 


5.40 



We see a useful feature of the N g = 4 results in Table 1: the Ck(t) are much 
more independent of the low-energy scale choice than in the N g — 3 case. This 
can be understood from Eq. (7). If we make a change in «&(£), then one can easily 
compute that: 



ACfc(t) 2Aa k (t)a k (t)a- 2 (t u ) 
C k (t) ~ ^ ~ 



a 2 (t v ) 



- 1 



(8) 



First, consider k = 1,2. Referring back to Fig. 1, we see that to a reasonable 
approximation Aai^ for N g = 3 and N g = 4 are very similar in magnitude in going 
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from mz to 1 TeV. Also from Fig. 1, we see that for k = 1,2 the denominator 
of Eq. (8) is very substantially smaller, while (Xk( m z) / a (tu) i s much larger, for 
N g = 3 compared to N g = 4. In combination, these two effects lead to negligible 
change in Ci ; 2 in going from mz to 1 TeV when N g = 4, compared to a modest 
change for N g = 3. For k — 3, the denominator is somewhat larger for N g = 3 than 
N g = 4, but the numerator is very much smaller for N g = 4, due not only to the 
decrease of a(tu)~~ 2 by a factor of ~ 4, but also because A«3 (in going from mz to 
1 TeV) is very much smaller for N g — 4 than for N g — 3 (see Fig. 1). A substantial 
change in C3 for iV^ = 3 is turned into a small change in the N g = 4 case. This 
is amusing both in its own right, but is particularly important for squark masses 
in that the small scale-sensitivity of C3 means the exact low-energy scale down 
to which we evolve will have a weak impact upon the squark masses. In the case 
of N g = 3, when determining the running mass of a heavy squark, it is crucial 
to evaluate C3 at the squark mass scale; for a heavy squark Cs(t~) can be very 
substantially smaller than the m^-scale value as seen in Table 1. In contrast, even 
for N g = 3, Cip are much less sensitive to a change in scale. 

Table 2: We tabulate the fik and Di values for squarks and sleptons. Also given are the mz- 
scale Ci values for both N g = 3 and N g — 4. For N g — 4, we employ m t ' = my — 100 GeV, 
m T i = m v i = 50 GeV, and tan/? = 1.5. We use the notation xw = sin 2 8w- 





fik 


A 


Q(Ng=4) 


Ci(N g = 3) 


i \ k 


1 


2 


3 








V 


1 

4 


1 





1 

2 


0.41 


0.52 




1 








— xw 


0.13 


0.15 


h 


1 

4 


1 





— \ + xw 


0.41 


0.52 


dR 


1 

9 





1 


~\ x w 


4.43 


7.32 


UR 


4 
9 





1 


\xw 


4.48 


7.37 


dh 


1 

36 


1 


1 


—\ + \xw 


4.80 


7.79 




1 

36 


1 


1 


2 — 3 Xw 


4.80 


7.79 



Let us now return to the squark and slepton masses. The primary sensitivity 
of these masses to N g is through the and the N g = 4 constraints on mo and 
mi/2- To see this let us recall that, in terms of the G\, one finds 

m\ = ml + Ciml/ 2 + Dim\ cos 2(3 , where Cj = ^ f ik C k . (9) 

k 

The fik and Di are tabulated in Table 2, as are t^-scale values for the C{, for 
both N g — 3 and N g = 4; in evaluating the CVs we have employed full two-loop 
evolution for all of the c^'s in determining the Cjfc's, as noted earlier. Also noted 
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above is the fact that, for N g = 4, employing the m^-scale values for the C% should 
be a reasonably good approximation, even in the case of heavy squarks, due to the 
modest sensitivity of the coefficients to the energy scale. Of course, we must 
keep in mind that the numbers given above have presumed the universal value 
Mi = m 1 / 2 at Mjj. Clearly, if the Mi do not have a common Mjj-scele value, one 

would have to redo all the computations. 

What about the value of (5 — tan -1 V2/ 'v{l Whether or not the slepton and 
squark masses are measurably sensitive to f3 depends upon the relative size of the 
itiq and C^m 2 , 2 terms compared to m? z cos 2/3. Certainly it would be nice to have 

sensitivity, since then it could be determined if tan/3 falls in the 1 ^ tan/3 ^ 
3 domain required if N g = 4. There are two mass-squared differences that, in 
principle, allow a direct determination of cos 2/3 independent of whether the mo's 
are universal and independent of any knowledge of the numerical values of the Cf. 

m J L - m l = m \ - m l L = - cos2/3m|(l - x w ) , (10) 

where we employ the notation xyy = sin 2 ^. (Note that since tan/3 > 1, cos 2/3 
is always negative and these mass differences are positive.) For such differences to 
be sensitive to cos 2/3 at just the la level, one must be able to measure a typical 
m to accuracy Am/m < m^/2m 2 . For masses of 100, 250, 500, 1000 GeV, this 
means Am/in < 32%, 5%, 1.3%, 0.32%; clearly the last two accuracies are at a very 
difficult level even for an e + e~ or collider of adequate energy. Determination 

of cos 2/3 at the 3a level of accuracy would require factor of three smaller errors than 
those listed above. To be more precise regarding the possibilities for measuring the 
mass differences in Eq. (10), we outline some relevant issues and possible techniques 
in Appendix B. 

The conclusion from Appendix B is that it will be difficult to measure mj and 
m~ to much better than about 5% accuracy. (Since the mass scales for and 
m~ L are even larger than for the sleptons, and errors in the determination of their 
masses will be larger also (see Appendix B), it seems clear that it is the slepton 
sector upon which we should focus.) Given the 5% accuracy estimate for mj and 

m~ determination and comparing to the criteria of the previous paragraph, we 
conclude that a direct determination of cos 2/3 will be problematical, especially if 
the slepton mass scales are as large as predicted for N g = 4. The ability to de- 
termine tan/3 could be even worse if tan/3 is large, since cos 2/3 varies slowly with 
j3 once j3 approaches n/2. Of course, for N g — 4 it is true that tan/3 is in the 
1 ^ tan/3 ^ 3 range where our ability to extract tan/3 from a cos 2/3 measurement 
would be maximal. Nonetheless, for the moment we must conclude that experi- 
mental determination of cos 2/3 and thence tan /3 will be difficult unless the slepton 
masses are well below 200 GeV, as is possible only if N g — 3. Of course, further 
study is undoubtedly warranted and could reverse this conclusion. Fortunately, 
the uncertainty in cos 2/3 is not the limiting factor in our ability to test the other 
crucial mass sum rules and relations discussed below. Other experimental and 
theoretical uncertainties are much more important in determining the limitations. 
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Sum rules that do not mix slept on/ squark and gaugino masses 



We can imagine two basic outcomes for our experimental attempts to determine 
cos 2(3: i) we are successful in getting at least a rough determination of m? z cos 2(3; 
or ii) our errors are too large to be sensitive to this kind of term. In either case, 
one would proceed next to a series of additional sum rules and mass relations that 
would provide tests of universality and the general unification/RGE framework 
and/or tests for N g = 4 vs. N g = 3. For the moment, let us assume that both 
squarks and sleptons are observable, if not at an e + e~ collider then perhaps at a 
collider of very high energy' 311 If not, one must proceed quite differently; we 
will return to this alternate case in Sec. 6c. 

To check mo universality independently of N g , one must consider mass combi- 
nations that are independent of the Cj. The only simple example is 



2(m~ — m~ ) + (m~ — m~ ) + (m~ — m~ ) 

v U R d R > v d R d L J v l L l R > 

= —x w m 2 z cos2(3 = yX H /(l - x w )~ l (m~ - m~ ) 
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This relation would be violated if all the mo's appearing in Eq. (9) are not the 
same. Because of the large number of mass differences appearing above, to detect 
a violation of universality in just one of these differences at the Amo/mo = f level 
could require that the squark and slepton masses be measured with an accuracy 
that is & //4. From the discussion of Appendix B, we conclude that it will be 
difficult to obtain accuracy / that is much better than 15—20%. However, this level 
of accuracy is acceptable in the sense that we will find that tests of all other sum 
rules are more or less restricted to this same rough level of accuracy even if the mo 
were universal. This means that the limited accuracy of the test of mo universality 
will not dominate our ability to check various mass relations. Let us presume that 
universality is shown to be satisfied within the accuracy of measurement achievable. 
Then we can proceed to extract mo and test for N g = 3 vs. N g = 4 to about this 
same level of accuracy. 

There are various mass combinations that can be used to extract a universal 
mQ. However, we should note that the four quantities mo, Ci, C2? and m| cos 2(3 
cannot be extracted using just the three slepton masses — we must have some 
squark-sector measurements. One simple combination is 

ml = m \ R - 3(m| fl - m|j + 4x w (l - x w y X [ml - m|) . (12) 
Note that the C3 terms that cancel between m~ and m~ will be evaluated at very 

UR d R J 

similar mass scales, so that the cancellation should be quite precise. The C\ terms 
that cancel between m~ and — 3(m~ — m~ ) could in principle be evaluated at 

i R y U R d R J F F 
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somewhat different mass scales. As discussed previously, this would not matter in 
the N g = 4 case, but would lead to a small uncertainty for N g — 3; however, this 
uncertainty is much less than that from simple experimental errors in determining 
m~ — m~ . Indeed, if mo is small, then experimental errors in evaluating the 

UR d R 

right-hand side above could become a severe problem. We will not dwell further 
on this issue here. 

The primary sensitivity to N g = 4 vs. N g = 3 derives from the fact that the 
gauge coupling at unification, a(tjj), is approximately twice as large when N g = 4 
as it is when N g = 3. This leads to the possibly large difference between the N g = 4 
and N g = 3 values of C3, depending upon the energy scale of evaluation, illustrated 
in Table 1. However, it is not straightforward to exploit this sensitivity through 
mass relations involving squarks and sleptons alone, that is relations designed to 
eliminate direct reference to m^2- This is already apparent from the modest 

changes in the Ck and C{ (see Tables 1 and 2) in going from 3 to 4 generations. 
The very best that one can do is to consider ratios such as C3/C1. As seen from 
Table 1, if squarks are at a low mass scale (of order mz) then this ratio takes the 
values 33.5 for N g = 4 as compared to 48.3 for N g = 3; but if squarks have masses 
of order 1 TeV, then we compute C3/C1 = 36.5, not very different from the 33.5 
value for N g = 4. Experimental extraction of C\ and C3 can be done in a number 
of ways. An example is: 

Cim\ /2 =3(m| fl - m 2 ~J - 3x^(1 - %)^(m~ - m| ) 

2 ' ( 13 ) 

C 3 m 1/2 ={m~ L - m~) + -(m~ R - m~J , 

where the C\m^, 2 result was already used in Eq. (12). Accuracy of about 10% or 
better would be very desirable for the experimental determination of Cim^, 2 and 



^3 m i/2 ^ we are ^° nave a good chance of distinguishing N g = 4 from N g — 3, but 
would require determination of the m~ — m-x and m~ — mj mass differences 

to roughly 3%, a challenging task even at an e + e~ collider, as we have already 
discussed. 

Some mass sum rules are quite insensitive to N g = 3 vs. N g = 4. Such mass 
relations can be used to test the general RGE approach, somewhat independently of 
the generation issue. To illustrate, consider the sample relation (3.19) in Ref. [28]: 



2 2 ( C ~2 ~ 4^1) 

z, — ?77,~ = — 

eR (C 3 - |Ci) 



( m ~ — m~ ) 

V d R e R J 



+ cos 2[3m z 



(14) 

3^ (C 3 - |d) v 2" 9 ^ 



2 (C 2 ~ A 9 , 

"0^77^ 877T _ (o ~ ^W 7 ) 



* We have corrected a sign error in this equation. 
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We note that cos 2(3 < (for tan/5 > 1) so that the sign of the second term is 
opposite to the sign of the factor in brackets. Adopting our results from Table 1 
for the Cfc, we find at scale mz'. 

ml -ml = (0.052 or 0.062) (m 2 ~ - ml )- [(19.3 GeV) 2 or (19.6 GeV) 2 ] cos 2(3 , 

&R dR "R 

(15) 

for N g = 3 or N g = 4, respectively. Since this relation assumes universality among 
the mo, the relatively minor changes in the coefficients for N g = 4 compared to 
N g = 3 implies greater sensitivity to possible universality violation among the mo 
than to N g . We note that for accurate examination of this sum rule, an accurate 
value for C3 is absolutely necessary. As noted above, in computing C3 it is critical 
when N g = 4 to include the two-loop contributions to the running of 0:3 given that 
the one-loop beta function is anomalously small, being proportional to 9 — 2N g . 

Third-generation mass sum rules (4.5) and (4.10) of Ref. [28], which are in- 
dependent of the mo universality assumption, are unaffected by a fourth family. 

However, since they require that the b mass matrix is diagonal to a good ap- 
proximation (which is generally true for the third family since m^ is small), the 

analogous sum rules do not apply for the fourth family due to the large value of 
my. 

A generalization of the sum rule (4.6) of Ref. [28] for the splitting ml — ml 

El 12 

can be found in the four-family case. It is not especially useful, however, since it 
relates ml, — ml, to m~ — ml via Af and Ay dependent terms; thus, we do not 

t l t 2 b 1 b 2 

display it here. 



Sum rules that mix squark/slepton masses and gaugino masses 



So far, we discussed only mass relations and sum rules that referred either 
to squark/slepton masses only, or to gaugino masses only. However, it should 
be clear from all our previous discussions that the real sensitivity to N g = 4 vs. 
N g = 3 lies in the dramatic shift of slepton/squark masses at given m-g. This was 
already discussed in Sec. 6a, and is further evident from the relations of Eq. (13), 
once the substitution ro^ = m~a(tjj) / a^(m~) is made and it is recalled that 
a(i[/)|jv ff =4 = 0.092 compared to ct(tu)\N g =3 = 0.043. However, because of inherent 
inaccuracies in experimentally measuring m~ it is better to reference M\{tz) or 
M2(tz), see Eq. (2), which can be quite well determined once the experimentally 
more accessible masses m^o and m~+ are measured. A particularly simple example 

is provided by m~ — ml . From Table 2 we find that (using m^-scale C{ values - 

for 1 TeV scale values the difference in coefficients is even greater) and the relations 
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ofEq. (2), 



- ) ■ ; - 1 ■ 1 2 

m~ -m~ =5m + (4.3 or 7.17)m 1 /2 + -x^/m^ cos 2/? 

=<5mo + (31.5 or 11.8)Aff(tz) + -x w m 2 z cos2(3 , 



for = 4 or A^j = 3, respectively. Here, we have included a possible violation 
of mo universality, <5mQ = rn^dji) — TOq(^r). Experimental uncertainties in deter- 
mining the coefficient, Cm 2 i °f M^tz) will be dominated by experimental errors 
in determining m~ — m~ and 5mk. The latter error, as discussed with regard 

6 d R i R ' 6 

to Eq. (11), is expected to be several times as large as the former. Very crudely, 
for experimental uncertainties in the squark/slepton mass measurements of size 
Afh/fh ~ /, one expects an experimental uncertainty in the M| coefficient of 
order ACm 2 ~ 3 — 5 x /to 2 /M|. We have already seen that / is very unlikely 
to be better than / ~ 0.1. Taking m ~ 1 TeV and M 2 ~ 200 GeV implies a 
ACm 2 ~ 7.5 — 12.5. This is smaller than the predicted N g = 4 vs. N g = 3 differ- 
ence, suggesting we would, in fact, be sensitive to N g by measuring masses of the 
gauginos and the squarks and sleptons of just the first two generations. 



6c. Non-Squark Spectrum Tests for Universality and Four Generations 

As repeatedly noted, mass determinations will be most easily performed at an 
e + e~ collider of appropriate energy. Regardless of the value of N g , chargino pair 
production is very likely to be within reach of a y/s = 500 GeV NLC. In the case of 
N g = 3, slepton pair production is also very likely to be possible for yHs = 500 GeV 
if boundary conditions of the no-scale or dilaton type are appropriate. For N g = 4, 
however, these same boundary conditions predict much larger slepton masses, and 
a higher energy machine would probably be required. However, for either N g — 3 
or N g = 4, the squarks are much heavier still, and it is quite possible that an e + e~ 
or collider would be built that could produce the neutralinos, charginos 

and sleptons, but not the gluino and squarks. Thus, it is desirable to consider 
what kinds of tests of universality, the RGE framework, and N g can be performed 
without reference to either the gluino or the squarks. In particular, we wish to 
determine if there is a test for N g = 3 vs. N g = 4 that can be performed using 

only the x?, X\ masses and the masses of the first or second family sleptons, which 
are those most likely to be both easily accessible and most precisely predictable 
{i.e. without reference to A parameters and L — R mixing). So, let us imagine 
that the M\ for i — 1,2 have been determined with reasonable accuracy from the 
observed values of m~o and m~+, following the procedure described earlier. If we 

Al Al 

then assume a value of N g , m^ 2 will also be computable from Eq. (2). Finally, 
the Ci can also be computed for the assumed value of N g . 
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As mentioned earlier, the £l and £r masses are very sensitive to N g when 
expressed in terms of the parameter £q- 111 what follows, we keep track of a possible 
difference between the £o values in the L and R sectors, i.e. of possible universality 
violation. We have, from Table 2, 



2 I>£l2 2 
771 ' 1 



7 ~[fo 1 m i/2 + (°- 41 or 0.52)m2 /2 + (-i + x w )m% cos 2/3 , 

(17) 

m f M£o¥ m i/2 + (°- 13 or 0.15)m^ /2 - x^m| cos 2(3 , 



for iV,, = 4 or Af„ = 3, respectively. We note that the C{ for mj and mj arc 
quite different, but not terribly sensitive to N g — 3 vs. 7V S = 4. However, by 
converting to the value of, for example, M\{tz) using Eq. (2), a large difference 
between the N g = 4 and N g = 3 phenomenology emerges. In the approximation 
that 777~o ~ M\{tz) we have 
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7771 
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>(29 or 7)[^] 2 + (12 or 3) + {-\ + cos 2(3 , 



9 

777~o 
X? 



(18) 



(29 or 7){^} 2 + (4 or 1) - cos 2/5 , 



777' 



777~ 
X? 



for N g = 4 or N g = 3, respectively. 

We plot our precise numerical results for rj and function of £0 m 

Fig. 16. We have taken m t = 165 GeV and tan/3 = 1.5. In the case of N g = 4, 
we choose = my = 100 GeV and 777 r / = m v t = 50 GeV. (All masses are 
running masses in this discussion.) The graph is given for the specific value of 
777~o = 117 GeV. For N g = 4 (3) this corresponds to 777 1/2 ~ 600 GeV (~ 300 GeV), 

respectively. This kind of plot has a number of important advantages. First, there 
is no sensitivity to the A parameter for the first and second generation masses 
being considered. Second, tan/3 enters into the location of the £0 = intercepts, 

but only weakly if m~o ~ mz- Third, by measuring the and £r masses in units 

of fn^o, the plot is quite independent of the actual value of 777^2 (or, equivalently, 

777-0 ), as seen in the approximate relations, Eq. (18). 

Fig. 16 exhibits a rather big difference for the functional dependence of rj 
and rj on £q> depending on N g . Suppose that experiment yields rj = 4.5. For 
N g = 3 this implies £q ~ 1.7 and, if £q = £q*, i.e. the £0 are universal, rj ~ 4.7; 
for N g = 4 we find £q ~ 0.8 and rj ~ 5.3. Since m~o ~ 117 GeV in this case, 
this corresponds to mj^(N g = 4) — mj^(N g = 3) ~ 70 GeV compared to perhaps 
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10 — 20 GeV uncertainty in the £l mass measurement. Of course, to measure mj^ 
and mj in the mass range above 400 GeV, required in this example, demands a 
very high energy e + e~ or fi + {i~ collider. 

As one moves to lower £o values, the discriminating power of this procedure 
slowly increases. Suppose the experimentally measured value is rj =3. If N g = 3 

this would imply = 1.1, and ^ = would then imply rj = 3.4. In contrast, 

if N g = 4, rj = 3 corresponds to £q = 0.45, yielding = 4.1, or mj^(N g = 

4) — mj L (Ng = 3) ~ 82 GeV, a bit larger than our previous result. 

How sensitive are we to a breaking of universality and the uncertainty in tan (31 
Consider first the possibility that £q and £q are not the same. Returning to our first 
example, and neglecting effects of the Tr(Ym 2 ) evolution terms, a measurement of 
rj = 5.3 in agreement with the N g — 4 prediction without universality violation 

could be reinterpreted as a prediction of N g = 3 provided £q ~ 1.9, as compared to 
the value £q = 1.7 required for = 4.5 for N g = 3. This translates to a breaking 
of universality in the amount: 

to 'to - °- 2 11%. (19) 



Thus, ~ 10% universality violation could correct for the N g = 4 vs. N g = 3 dif- 
ference. This is an intrinsic ambiguity when only slepton, chargino and neutralino 
masses are available. 

In the above, we have purposely chosen situations where consistent solutions 
for both N g = 3 and N g = 4 are possible. It could easily happen that consistent 

choices for £q and £q are simply not possible for the observed values of r~ and/or 

rj . In particular, if rj ^ 2.0 or rj ^ 3.5 then N g — 4 is excluded, whereas N g — 3 

models would be possible. In fact, if N g = 3 and £o is in the preferred £o ^ 1 range, 
then rj < 3 and N g = 4 would be immediately excluded. The lower bounds on 

rj^ and rj for N g — 4 are, of course, equivalent to lower bounds on mj^mj 

of 410, 234 GeV for this particular m~o value. And, these lower bounds can only 

be reached if one is willing to accept the no-scale itlq = boundary conditions. 
In another extreme, it could happen that one measures r-j ~ rj ^7 — 8. Such 

Zr Zl 

large slepton masses relative to m~o require very large £o values if N g = 3, well 
beyond any acceptable model range. In contrast, if N g — 4 the required £o values 
are perfectly reasonable. 

We note that sensitivity of these results to cos 2(3 is really quite limited. For 
m~o ^ mz, and taking both xw and —\ + xyy to be of magnitude ~ 0.25, the 
full range of cos 2/3 between and —1 corresponds to Ar ~ 0.125/r x mz/m~o ^ 
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0.125/r. Since the r values are typically at least of order 1, most probably much 
larger, we see that the cos 2/5 uncertainty is, at worst, roughly the same as the 
experimental error, and is not significant compared to the predicted r differences 
or possible universality questions. 

Finally, we note that even if we have an ambiguity in the decision between 
non-universality and N g = 3 vs. N g = 4, in general we will have a very rough 
determination of the general size of £o- m the above example, we would at least 
know that the £o values were ^ 0.8 if N g — 4 or ^ 1.7 if N g = 3. This would allow 
us to set the overall boundary condition picture. 



7. Experiment and a Fourth Family 

In this section, unless otherwise indicated, all masses are pole masses. Where 
appropriate, running masses will be explicitly written in the m(m) notation. Very 
roughly, pole masses are higher than running m(m) masses by about 5 or 6%, 
based on the dominant correction m(pole) = m(m)[l + ^a s (m)]. 

In all our discussions up to this point, we have adopted the conventional 
interpretation that the CDF and DO events derive from t — > bW decays, with 
mt(pole) ~ 175 GeV. Alternative possibilities in which the events at CDF and 
DO are reinterpreted as my ?Z 175 GeV events with t' — > bW, while the t is not 
currently observed, are not consistent within the perturbative MSSM GUT frame- 
work. In the first such option, the t is not observed because it is too heavy, 
mt(mt) 3> 165 GeV. The difficulty with this option is that for such large values 
of mt the maximum value of mt'(mt') allowed by perturbativity is much below 
165 GeV and the t't' events would not have the same characteristics (e.g. recon- 
structed mass and cross section) as observed by CDF and DO. In the second such 
option (emphasized in Ref. [14]), the t is much lighter than the t', mt ~ m^, but 

is not detected because its decays are dominated by t — > tiVi- However, in the 
MSSM GUT context, the latter decay does not occur if the soft-SUSY-breaking 
parameters are approximately universal, since (as we have seen) the t\ is always 
much heavier than the t. Further discussion of this latter case will appear in Sec. 
7a. Certainly, the phenomenology of the t' and b' at the Tevatron that we present 
depends crucially on the mt ~ 175 GeV > mv,mt' hierarchy, not to mention the 
relative size of my vs. my , and on the mixing pattern between the generations. 

In Sec. 7a we will discuss direct limits on the fourth family fermions coming 
from collider experiments, especially the CDF and DO experiments at the Tevatron, 
assuming that mtipole) ~ 175 GeV. We find that a relatively small portion of 
my , m t > parameter space survives all such constraints for the more natural choices 
of CKM mixing between the fourth generation and lower generations. We also 
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describe in more detail why the mt ~ m\y scenario is not viable for universal soft- 
SUSY-breaking boundary conditions. In Sec. 7b we give the constraints on mp, my 
arising from consistency with precision electroweak data. Not surprisingly, these 
tend to rule out large mt> compared to my. Finally, in Sec. 7c we speculate on 
deviations in single-jet inclusive and di-jet spectra that could arise as a result of a 
full four-family content at relatively low mass scales. 



la. Direct Collider Searches 

As we have demonstrated, the requirement that all Yukawa couplings remain 
perturbative places rather strict upper bounds on the t', b' , v' and t' masses, and 
the possibility of an MSSM fourth generation may well be settled within the next 
few years by searches (i) at the Tevatron, for an appropriate set of new signatures 
(i.e. besides t — > b + W), and (ii) at LEP-II, via detection of direct r' + f' and/or 
v' + v' production. 

LEP has set the firmest direct limit, namely mj ^ mz/2 for / = b' ,t' , r' ,v' . We 
have also seen that the value of mj ~ 175 GeV, as determined by the recent CDF 
and DO results) 201 is such as to imply mj > my , me for perturbative consistency. 

Let us begin by noting the critical difference between the mt> < my and mp > 
my cases, assuming that mt > mt>,my. If mp < my, then the t' — > Wb decay 
proceeds at tree-level at a rate determined by the CKM entry Vpy One-loop 

processes, such t' — > eg, t' — > c~y, t' — > cZ^ and so forth occur at rates determined 
by Vt' s V cs and VpyV c y times a loop integration factor. It would be extremely 
unnatural for the product of a loop integration factor times a two-generation- 
skipping CKM matrix element to exceed the one-generation-skipping CKM Vpb- 
Thus, it is almost certain that t' — > Wb decays would be completely dominant if 
mp < my. This situation should be contrasted with the reverse case of mp > 
my, for which the interesting issue is how the b' decays. The decay b' — > Wc 
occurs at tree- level with rate determined by V c y. The one-loop rates for b' — > bg, 

b' — > 67, b' — > bZ^ and so forth are determined by V t yV t b and V t iyV t >b times loop 
integration factors. It is actually rather likely that the two-generation-skipping 
CKM element V c y is smaller than the one-generation-skipping Vty or Vyb times 
the loop integration factor. In this case, the flavor changing decays b' — > bX would 
dominate over the two-generation skipping b' — > Wc decay. The different decay 
structure for wif < my as compared to m? > my has crucial phenomenological 
implications. 

To repeat, if m t > < my it is almost certain that the t' will decay in a top- 
quark-like manner to Wb (or W*b, depending on mp). Events of this type are 
clearly ruled out by early CDF and DO 'top-quark' limits for all mp < my regions 
allowed by the four-generation perturbativity constraints, see Fig. 2. Thus, we 
immediately eliminate roughly half of the perturbatively allowed regions. We now 
focus on experimental restrictions that arise if m t > > my. 

b'b' Production at the Tevatron 
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If rrif > my, then in b'b' production at the Tevatron we must consider the two 
more or less competitive b' decay scenarios outlined above. If Vy c (or, much less 
likely, Vy u ) is sufficiently large that b' — > cW (or b' — > uW) decays are dominant, 

[32] 

the published data already rule out my < 85 GeV at 95% CL. (Presumably this 
old result will soon be updated using the full CDF and DO data sets.) But, if the 
b' is essentially unmixed with light generations, it will decay via flavor- changing 
neutral current (FCNC) channels: b' — > by, bg or bZ* for my < mz, with b' — > bZ 
becoming dominant for my > mz + my 

Present data would appear to rule out a significant rate for b' — > bZ decays 
by virtue of there not being room for an excess number of Z's beyond that pre- 
dicted on the basis of the standard qq annihilation mechanism. Taking a b' mass 
of 100 GeV (roughly the maximum possible if my > my), the cross section (with- 
out inclusion of any K factor) at yfs = 1.8 TeV (Fermilab Tevatron) is about 
60 pb. For an integrated luminosity of 67 pb -1 this implies roughly 4000 b'b' 
events. If we imagine looking for Z's from the decay of either (or both) 6"s, then 
this number is to be multiplied by [1 — BR(b' — > non— Z) 2 \. Given the small phase 
space, my — mz — m^ ~ 5 GeV and the low value of mp ~ 100 GeV required by 
our perturbative constraints (which impacts the one-loop Penguin diagram cal- 
culations), we expect' 71 BR{b' — > non— Z) to be significant but probably ^ 0.7, 
implying > 2000 events in which at least one of the 6"s decays to a Z. Using 
BR(Z — > £ + £~) = 0.067, and an overall two-lepton cuts/acceptance/detection 
efficiency of about e = 0.5, we find ^ 80 events where at least one Z appears 
in the b'b' final state and has a clearly detectable leptonic decay. (If we take 
BR{b' — > Zb) = 1, and adjust for luminosity, this number is consistent with 
the pre-pT(Z)-cut numbers of Agrawal, Ellis and Hou, Ref. [7], for the sum of 
£ + £~ + jj + X and A£ + X event numbers at my = 100 GeV.) In addition, Fig. 2 
shows that if m t > is to be larger than my and my ~ 100 GeV, then m t ' can only be 
just larger than my. In this case, the t't' production cross section would be more 
or less the same as the b'b' cross section, thereby almost doubling the 80 events, 
assuming that t' — > b'W* decay is dominant. Although a specific analysis has not 
been presented for the latest CDF or DO data with regard to this point, we expect 
that such a large number of additional relatively clean Z events can be ruled out. 

For an integrated luminosity of L = 150 pb -1 , it would seem that one could virtu- 
ally eliminate the my > mz + m^ case down to very near the threshold. However, 
if b' — > bh° is also kinematically allowed, the number of extra Z events would be 
much smaller and my > mz + m^ would not be excluded. We shall return to this 
issue shortly. 

Thus, consistency with experiment probably requires that my < mz + m^, 
and that b' decays are dominated by the FCNC channels outlined above. These 
'non-standard' decays would yield distinctly different signatures than the lepton- 
plus-jet signatures of charged current decays that are already mostly excluded, 
and new search strategies are required. Perhaps the cleanest discovery channel 
makes use of b'b' production followed by double b' — > 67 decay. As an example 
consider my = 90 GeV(95 GeV). The Tevatron cross section (without including 
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any QCD K factor enhancement) for b'b' production is ~ 100 pb(76 pb), corre- 
sponding to ~ 6700(5100) events for L = 67 pb -1 . Assuming BR(b' — > try) ~ 0.1 
(typical of the low rrif results given in Hou and Stewart, see Ref. [7]), we obtain 
about 67(51) 6677 events. In order for these events to be sufficiently free of back- 
ground, single b tagging is probably necessary. If we assume that the efficiency 
for tagging of at least one b (after appropriate kinematical cuts) is approximately 

50% (as found in the most recent CDF top-quark analysis' 201 ) and that the overall 
cuts/acceptance/detection efficiency (after requiring minimum transverse momen- 
tum, say pt > 15 GeV, and central rapidity for the photons and 6's) is about 35%, 
we see that a clean sample of 11(8) or so events should be present in the CDF 
data. Of course, t't' production will also contribute to the 6677X signal assuming 
t' — > b'W* is the dominant decay of the t' . (For mt> near my the jets or leptons 
from the W* will be quite soft, and the signal rather similar to that for direct 
b'b' production). Events of the 27 + 2j type with at least one tagged b quark 
are currently under investigation by the CDF collaboration) 331 and we expect that 
CDF will shortly be able to severely restrict the allowed four-generation scenarios. 
Since event rates rise rapidly with decreasing my, it is hard to imagine that a b' 
with mass much below 80 — 90 GeV would escape detection. Referring to Fig. 2, 
this means that the my (my) < 75 — 85 GeV region for tan/5 = 1.5 and the entire 
tan/3 = 2.2 region are on the verge of elimination if no signal is reported. An 
integrated Tevatron luminosity of L = 150 pb -1 would surely allow one to use the 
6677 final state to eliminate all my ^ m^ + mz, extending up to my values such 
that the b' — > bZ mode really becomes significant enough to suppress the b' — > 67 
branching ratio. At this point, as described above, one looks for extra Z's, which 
also provide an excellent signal. 

Of course, if a reasonable K factor is included in the b'b' cross section (K = 1.3 

to 1.5), the event rates for the 6677 and/or extra Z events will increase above 
those given above. This leaves even less room for the four-family model to avoid 
detection. 

The only possible escape from the b' — > 67 and b' — > bZ decay mode constraints 
would be if b' — > bhP decay is kinematically allowed. If allowed, it will certainly 
dominate all other modes when b' — > bZ is kinematically forbidden. Even when 
the b' — > bZ channel is open, b' — > bh° will at the very least be competitive 
and suppress the b' — > bZ branching ratio — indeed, the bhP channel could still 
be the dominant mode.' 341 However, as seen in Fig. 6 (after converting running 
masses to pole masses), this escape is not available for the my > my portion of 
parameter space in the case of the dilaton boundary condition scenario, since there 
my < m^o . Only in the extreme lower right-hand corner where the running masses 
obey m/jo ~ my ~ 100 GeV and my — 50 GeV, is the b' — > bhP decay kinematically 
allowed, and this corner is almost certainly ruled out (as discussed in an earlier 
paragraph) since m t < < my. However, for scenarios with larger (but not too large) 
mo, the hP can be light enough to make the bh° channel kinematically accessible. 

In Fig. 11 we plotted as a function of £4, for various fixed £0 values, the 
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minimum m~(m~) value that is allowed by the r[ constraints, taking my (my) — 
mt>(mt>) = 100 GeV (corresponding to pole masses of about 105 GeV). We can 
now use this minimum m~(m~) value to fix all other model parameters, thereby 

determining the sparticle spectrum, in particular m^a . Since the m? Hi contribu- 
tions to the 777? = m? H . + jj? mass terms in the scalar potential will tend to be 
smallest when m~ is small (and \x also tends to scale with 777~), these minimum 
m~ values give m^o values that will be very close to the minimum allowed for the 

fourth-generation fermion masses chosen. To determine if b' — > hPb is kinematically 
allowed, we compute 7775/ (pole) — m^o — m&. Plots of m^o and this mass difference 
also appeared in Fig. 11. From the my (pole) — m^o — curves, we see that for 
middle-of-the-road values of £0 the b' — > bhP decay can be kinematically allowed. 
When the b' — > bhP decay channel is present, thereby suppressing the b' — > bZ 
channel branching ratio, a b' with pole mass somewhat above 105 GeV could have 
escaped observation. 



t't' Production at the Tevatron 

Let us now turn to a more thorough examination of t't' events. For the my > 
100 GeV mass region, there are two competing decay modes for the t'\ t' — > b'W^ 
and t' bW . The b'W decay will certainly dominate over the bW mode when 
the former W is real, since the bW channel is suppressed by the intra-generational 
mixing factor Vyb- However, when my — my < mw, the two-body bW mode could 
be competitive with the three-body b'W* decay. In Fig. 17 we plot the value of Vyb 
for which Y(t' — > b'W*) = T(t' — > bW) as a function of my for several 7775' values 
(we employ pole masses). For Vyb larger than the values plotted, the t' — > bW 
mode is dominant. We see that for b' (pole) masses larger than the rough lower 
bound of 80 GeV, discussed above, it is not at all impossible that Vyb could be 
large enough for t' — > bW to be the dominant mode. 

For my % 170 GeV (as required by our Yukawa perturbativity bound), such 
dominance would imply a large excess of events relative to those already present 
by virtue of it production and decay. As discussed in more detail below, the CDF 
and DO experimental results are now more or less consistent with the expected it 
rates for mt ~ 175 — 185 GeV. Taken at face value, these results then imply that 
t' — > bW cannot be the dominant t' decay mode. 

In fact, we show below that it is only possible to obtain sufficiently few CDF 
and DO events if the t' decay is to b' plus a highly virtual W*. A typical case 
that survives is my ~ 115 GeV and my ~ 80 GeV. From Fig. 17 we find that 
Vyb ^ 0.02 (a value roughly the same as V c b) is required for t' — > bW to be 
adequately suppressed for these mass choices. Such values for Vyb are certainly 
reasonable; in fact, one might expect considerably smaller values based on the fact 
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that Vfjc <C V us , but we must also consider implications for b' decay. We noted 
earlier that if b' — > cW were the dominant b' decay mode, then b'V events would 
have been detected. Consistency with experiment thus requires that the one-loop 
b' — > bX neutral current decays are dominant. This, in turn, requires (see Hou 
and Stewart, Ref. [7]) that Vt^/Vyc, V t y/Vy c ~ 10 2 - 10 3 . Combining with the 
V t 'b restriction above, we find that Vy c ^ 10~ 4 is required. As seen from Fig. 17, 
this limit would become stronger for larger values of my. Suppressions larger 
than those suggested by measured V c \, and V u \, values would be required for my 
values much above 100 GeV and thus become problematical. As discussed earlier, 
my ^100 GeV is also the range for which b' — > hPb decay dominance is probably 
required to avoid too many extra Z events from b' — > Zb decays in b'V events. 

In our previous work) 11 we performed a rough Monte-Carlo study of the number 
of t't' production and decay events expected in the dilepton-plus-jets and lepton- 
plus-jets channels, assuming that t' -> b'W^ is the primary t' decay. At the time, 
the number of additional events (beyond the predicted number from tt production) 
in these two channels was not inconsistent with the then-available CDF data for any 
of the larger m?, my choices within the perturbatively allowed region. However, 
the luminosity accumulated by both CDF and DO is now much larger [L ~ 67 pb _1 
vs. 19 pb _1 ). Further, CDF now employs new 6-tagging algorithms which roughly 
double their 6-tagging efficiency Consequently, the constraints on the model are 
now much stronger. Exactly how severe the constraints are depends on the precise 
predictions for event rates from normal ti production. The CDF collaboration 
states that their results are in good agreement with a top-quark Monte Carlo. As 
a cross check, we have also repeated the Monte Carlo study of our previous work, 
Ref. [1], after adjusting for the new luminosity and new 6-tagging procedures. 
Despite the approximate nature of our implementation of the CDF cuts and b- 
tagging procedures, agreement is good. For mt = 175 GeV, the uncut ti cross 
section is about 3.5 pb, and roughly 25 events pass the single-lepton plus b-t&g 
plus > 3-jets criteria, while ~ 3 events pass the dilepton cuts. CDF observes 37 
events in in the W+ > 3 jet channel in which at least one b is tagged, with an 
estimated background of 16, and 6 events in the di-lepton channel, with estimated 

background of 1.3. Thus, their number of signal events is 21 and 4.7 in these two 
respective channels, to be compared to 25 and 3. Clearly, there is little room for 
additional events passing their cuts from t't' production. 

The rates (after implementing the same cuts, etc.) in these two channels 
deriving from t't' production, followed by FCNC decays of the type b' — > b + jet(s) 
(in particular, no extra photons or leptons from b' decay are allowed) appear in 
Table 3. (In Table 3 we have included several (mt>,my) pairs that are excluded 
by the boundaries of Fig. 2 just to indicate the effects of different mass choices 
on the rates.) Restrictions on the perturbatively allowed regions of Fig. 2 are 
now significant. For tan/3 = 2.2, for which only values of my < 85 GeV (i.e. 



* The number of background events, 16, in the W+ > 3-jcts channel is estimated by scaling 
the tag background estimate of 22 events by the ratio 37/50 of events/tags. 
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Table 3: For L = 67 pb 1 , we tabulate the number of t't' events passing our approximations to 
the CDF top-quark discovery cuts in the single lepton+6-tag and dilepton modes. 





Lepton + 6-tag Mode 


Dilepton Mode 




my (GeV) 


my (GeV) 


mv (GeV) 


50 80 110 130 


50 80 110 130 


160 


20 16 3.6 1.5 


2.4 1.8 0.2 0.01 


130 


53 20 0.5 


4.2 1.1 0.002 - 


100 


58 3.6 - 


2.8 0.01 - 



my (my) ^ 80 GeV are perturbatively allowed, Table 3 shows t't' production will 
yield sufficiently few CDF single-lepton + 6-tag events only if my ^ 105 — 115 GeV, 
effectively eliminating the upper 50% of the allowed m t >(m t i) region. For m t i ^ 
105—115 GeV, Table 3 also shows that my can't be too light; roughly my ^ 70 GeV 

is required? The underlying reason for these restrictions is easily summarized. For 
wif ^ mti the only way the t't! events can evade being included in the CDF event 
sample is if m t > — my is sufficiently small that the W in t! — > Wb' is virtual and 
the jets and leptons from the two VF*'s are quite soft. For such cases, additional 
t't' events could be included in the event sample if CDF and DO could soften their 
cuts. However, backgrounds probably increase dramatically. If this is the case, 
then additional luminosity is the only way to improve sensitivity; for example, 
L = 150 pb -1 would yield about 8 CDF events (for current cuts and 6-tagging 
efficiency) if m t > = 100 and my = 80. This would probably constitute a detectable 
signal. 

Combined Restrictions 

Combining the t' restrictions with the existing and probable b' CDF data con- 
straints, and accounting for the limited range allowed for my by perturbativity (es- 
pecially for tan/? = 2.2) and the m t > > my requirement, we are clearly forced into 
the running mass domain 110 GeV ^ mt>(mt>) ^ my (my), 95 GeV ^ my (my) ^ 
75 GeV. For tan (5 — 1.5 this is the region that is as close to a fixed-point limit as 
the perturbativity constraints allow and is also the region for which the unification 
prediction for a s (mz) is smallest and, thus, in best agreement with data. We have 
also seen that it is a region for which Vfb and Vy c can both take on values that 
are reasonable extrapolations from those already measured in the lighter quark 
sectors, while at the same time being consistent with dominance of t' — > b'W* over 
(probably excluded) t' — > bW decays in t't' events and of b' — > Xb neutral current 
decays over (probably excluded) b' — > cW decays in b'b' events, provided m t » (my) 
lies in the upper (lower) portion of its allowed region. A tripling of the Tevatron's 



f Note that the t' — > b[xt decay, that could provide an escape from these restrictions, is not 
kincmatically allowed for any of the mass spectra scenarios of Figs. 13-15. 
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integrated luminosity to L ~ 150 pb would almost certainly either close the my 
window above, or lead to b' discovery 

The Hidden- Top Scenario 

To escape from the experimental constraints on t'i' production and decay 
discussed above would require that the 'top-quark' events not come from third- 
generation ti production. As described earlier, in one such scenario mj ~ mw 
and t — > t\x{ decays are dominant. Typical (pole) masses required for this sce- 
nario to be viable are those discussed in Ref. [14]: m t ~ m^, m t > ~ 170 GeV, 
my ~ 110 GeV, m~o ^ 25 GeV, and mj ^ 60 GeV. We have explicitly analyzed 
this scenario using our full SUSY GUT framework. We find that for m T > = 50 GeV 
(the choice of Ref. [14]) and m v i ^ 78 GeV, the above mt, mt> and my masses are 
consistent with gauge coupling unification and with perturbative evolution for the 
Yukawas up to the scale M\j . (Our criteria do not quite allow consistency with 
Yukawa perturbativity for the m v i = 80 GeV choice of Ref. [14].) 

However, we find that it is not possible to obtain m~o + mj < mt (where t\ 

is the lighter stop mass eigenstate) in the context of SUSY GUT scenarios with 
universal mo and A parameters at M\j . Since both m^o and scale with m~, the 
minimum value for the mass sum will occur at the smallest possible value of m~. 
As described earlier, the minimum m~ (or equivalently mi^) value is set by the 
t[ constraint. We find that the r{ constraint (for the mass scenario just outlined) 
implies an absolute lower bound on the M(y-scale parameter m^2 that is similar 
in nature to that illustrated in Fig. 10, except that lower values of m^^ (close to 
100 GeV) are allowed at high £o- This translates into m~ ^ 150 GeV. Analogously 
to Fig. 11, much higher values are required when £o = m o/ m i/2 ~ 1- Also as in 
Fig. 11, the minimum value of m\ji is achieved for a range of £4 = Ajm^n values 
centered about which is broad for large £o> but increasingly narrow for £0 ^ 1. 
However, even at large £o 5 arbitrariy large values of \^a\ do not yield consistent 
solutions for any as illustrated by the termination of the low-£o curves in 

Fig. 11. Too large a value for |^| (typically |d|/£o ~ 4 — 5 at large £0) results 
in either m 2 A0 < at the EWSB minimum, and/or m 2 < for one of the colored 
and/or charged 4th-family sparticles. 

We describe two typical cases: £0 ~ 1 an d Co ~ 5. For £0 ~ 1> m g ~ 400 GeV 

is required, with m~o ~ 0.125 — 0.13m~ implying m~o <: 50 GeV. For £q ^ 5, we 

ai y Ai 

find m~ 150 GeV and m~o ~ 0.14m~, implying m~o ^21 GeV. ( Both of these 

minimum values only apply for small £a/£o- As £a/£o increases, the minimum 
allowed value of m~ increases, but very slowly in the case of £0 — 5.) Although 
the latter m~o value is in an acceptable range for the scenario, is always quite 
large for £4 values that do not imply m 2 A0 < and/or color/U(l) breaking. This 
is because is always well above mt = 80 GeV at small £a (typically by three 
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hundred GeV or more for £o = 5). Thus, even though ti — £r mixing increases 
so that mj decreases as one moves to higher in the ^ < direction, mj 

only declines by some 50 GeV before \£a\ is large enough that m 2 A0 < and/or 
color/U(l) breaking arises. 

For allowed £4 values, what prevents a sufficiently small mass for the t\ is 
the fact that the r{ cannot be too light. The simplest way to avoid this lock-step 
arrangement is to break the universality assumption for the tuq parameters at Mjj, 
allowing much smaller itlq in the t sector than in the r' sector of the theory. We 
have not explored this option. 



7b. Four- Family Precision EW tests* 

In addition to direct collider searches for the presence of a possible fourth 
family, indirect effects arising from virtual quantum effects can be probed through 
precision measurements of many observables and can lead to indirect limits or 
bounds on the parameters of the theory. For example, the well-known p-parameter 
limits the magnitude of iso-spin breaking, and thus bounds the m\, — m\ and 
m\, — vr? v , splittings of a possible fourth family. On the other hand, although 

the LEP measurement of Z — > bb is a particularly sensitive probe to any new 
physics involving possible coupling to the top quark, the vertex contribution to 
Z — > bb arising from a fourth family is expected to be small due to CKM mixing 
suppression. We therefore neglect all vertex contributions in our analysis, focusing 
instead on vacuum polarization effects. 

Overall, it desirable to perform a global fit to the present experimental data 
using several different observables. The most comprehensive approach would be 
to perform the complete, one-loop calculation for each observable in question, and 
do a best fit of the new physical parameters of the model, (i.e. mt>, my, m T ', 
t?V, tan/5, /x, . . .). This rather ambitious approach has only recently been at- 
tempted for the MSSM , and will not be attempted here. A more modest approach 
inevitably involves making certain assumptions and approximations, i) As previ- 
ously explained, a fourth-family added to the MSSM contributes dominantly via 
'oblique', or vacuum polarization effects, ii) We assume that the contributions 
from the sparticles are small enough to be neglected. This is likely to be a good 
approximation for N g = 4 models. Indeed, we have seen that for N g = 4 the su- 
pergravity boundary conditions and evolution equations imply a fairly large lower 
bound on m~ (especially for the dilaton and no-scale boundary condition choices) 
which, in turn, implies that most of the sparticle masses are large. According to 
the decoupling theorem, loops involving massive SUSY particles contribute negli- 
gibly to vacuum polarization amplitudes. The particles most likely to violate this 
approximation are the sometimes- light r[ and u[, and even these are heavy for 



★ In this sub-section we will exclusively employ running masses. 
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larger m-x values, independent of model. Finally, iii) we assume that the h° of the 
Higgs sector is SM-like, implying that the SUSY Higgs sector is equivalent to a SM 
Higgs sector with a light Higgs boson. Given the four-family constraint tan/3 ^ 3, 
this is true for m^o ^ 200 GeV, as is essentially always the case for the models 
considered. 

[35] 

In the following analysis, we perform a global fit to the latest LEP+SLD data 
and we employ an extension of the Peskin/Takeuchi S,T,U formalism 1361 in order 
to find the 90%, 95% allowed regions of parameter space. In this approach, for 
a given 'reference' SM (i.e. with chosen mt,mn values), S,T,U = 0, and the 
(non-zero) best-fit values to experimental data are a measure of the oblique, or 
vacuum polarization, contributions that should be explained by 'new physics'. For 
our reference SM we adopt N g — 3 with m t (m t ) = 165 GeV and Higgs mass 
run — 100 GeV (typical of the m^o values found in our computations). We then 
treat the contributions from t', b', v', t' as new physics. 

At this point, there are several parameter choices and assumptions that affect 
our results, i) What values of a em (mz), a s {mz) are chosen for the reference SM? 
As was recently shown, new estimates of a em (mz) can shift the best-fit for the 
experimental value of the S'-parameter relative to S = for the reference SM 
by +10%. Similarly, the choice of a s (mz) is quite important. Secondly, ii) one 
should be careful when using the S, T, U formalism in the presence of new light 
physics, since in this case it has recently been shown that the standard S, T, U 

formalism is simply not an adequate parameterization.' 37 ' 381 In the original work 
of Peskin and Takeuchi, the S, T, U parameters were explicitly defined using the 
linear approximation in an expansion in g 2 /M|! 361 Although perfectly valid when 
considering new physics with scales much higher than mz, such as technicolor, 
this expansion is not accurate for new physics near or below mz- It is possible to 
define a set of parameters which do not rely on a g 2 /M| expansion. In a complete 
comparison to all available experimental observables, this results in a proliferation 

of new parameters that must be simultaneously fit? However, if one restricts the 
'global' fit to Z-pole plus m^/m^ measurements, then one can show that a re- 
definition of S — > S', T — > T', U — > U' is possible, where the primed variables can 
be computed in the presence of additional light physics without approximation. 
Further all the observables in question are linear functions of S', T', U'. Thus, the 
procedure is to restrict the fit to measurables which are defined at the Z-peak (along 
with the mw/mz measurement), employ the standard three-parameter S,T,U 
global fit to the experimental data (relative to a reference SM with S = T = U = 
by definition), and finally re-interpret the best fit S,T, U values as S',T', U' . 

In order to compare to predictions of a four-generation model, we calculate the 
exact (no expansion is performed) contributions to S', T', U' from loops involving 



f In a quadratic expansion, in general seven independent parameters arise ' , and in a fully 
exact treatment, eleven parameters must be specified 1391 
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fourth-generation particles for a given set of m<, ra<' , my : m u ' : m T > values; denote 
these contributions as S'^T'^U'^. Since the observables O are linear functions of 
S',T',U', we can then compute Oi(N g = 4) = O l (S' 4 ,T^,U' 4 ). The predicted 
values of Oi(N g = 4) are then compared to the best fit experimental values as 
obtained using ZFITTER (relative to the reference SM assumed) and an overall 

Ax 2 computed. This computation is performed for each point of interest in the 
fourth-generation parameter space. As noted earlier, we use a reference SM with 

int{mt) = 165 GeV and m# = 100 GeV and take a^^mz) = 129.08 [1T1 . We employ 

the latest values for the following measurements: mw/mz, T z, Ri, Rb, R& ^fbi 

A fb> FBi A T , A e from LEP and A° LR from SLD. [351 Finally, because of sensitivity 
of the reference SM prediction from ZFITTER for these observables to a s , we give 
results for the two values, as(mz) = (0.12,0.13). 

The solid and dashed lines in Fig. 18 summarize the final results of this global 
fit. There, we have adopted m T ' = m v i = 50 GeV, tan/9 = 1.5 and have allowed 
nT>t'{ m t') and mb'(mb>) to vary. The region in the m t ',mb' parameter space above 
the lower (upper) solid line is excluded at 90%(95%) C.L. if the reference SM has 
ot^{mz) = 0.13. The two dashed lines give the same C.L. boundaries in the case 
a z( m z) — 0.12; in this latter case, Fig. 18 shows that there is a region at small mf 
that is also excluded. As one might have naively anticipated, the region of large 
m? compared to has a poor C.L.; however, we also see that if both and 
m? are small, then predictions of the fourth-generation model for S',T',U' can 
again deviate significantly from the experimentally preferred values, depending 
upon a s (mz)- The region of modest m t '(m t ') ~ mb'(mb>) ~ 100 GeV, preferred on 
the basis of current and anticipated Tevatron results, lies well within even the 90% 
C.L. limits. 

7c. Influence on jet and di-jet spectra of a slowly running four-family a s 

One other indirect effect of the presence of a fourth family, along with its 
full complement of supersymmetric partners, all at relatively low energy scales, 
has recurred throughout our discussion: a s (Q) decreases much more slowly with 
increasing Q than in a standard three-family model, especially if the comparison 
is with the case where the masses of the superpartners of the three-family model 
are taken to be large. Thus, it is important to not forget that a s (Q) is, in fact, 
directly measurable through the single-jet and di-jet spectra of light quark jets. Of 
course, the exact scale Q at which a s is evaluated will be an important issue. In 
the MS scheme, higher order calculations suggest that an appropriate Q value for 
single-jet inclusive spectra and di-jet mass spectra is the subprocess center-of-mass 
energy, roughly given by 2E^ and rri2j, respectively. How large can the deviations 
be? To illustrate, we consider as a function of Q the fractional enhancement ratio 

MV g = 4)] 2 -[af^(iV g = 3)] 2 
Qs " [aS M (N g = 3)\ 2 ' 1 } 

where ctf M is the a s predicted in the three-family case with no superpartner effects, 
(i.e. with all superparticle masses taken to be large); results for a s (N g = 4) will 
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be illustrated for my {pole) = mf{pole) = 105 GeV (our preferred mass range) and 
^SUSY va l ues °f 200 GeV and 400 GeV. Such values are quite typical of the overall 
mass scales for the strongly-interacting sparticles (that contribute to a s evolution) 

in the model scenarios we have studied?" We employ a 2 s in Eq. (20) because this is 
what appears in the QCD subprocess cross sections. 

The results appear in Fig. 19. We see that for M^^y = 200 GeV, percentage 
deviations are, for example, about 25% at Q — 1 TeV. This would correspond to 
subprocess energies of 1 TeV for the standard MS scale choice mentioned above. 
The corresponding value for the (probably more reasonable) choice of M^^y = 
400 GeV is 12%. If there are only three families with superpartners at a light mass 
scale, the values of R aa , Eq. (20) ( with a s {N g = 3) in place of a s {N g = 4)) are 
very much smaller in all cases. Thus, it is not entirely impossible that preliminary 
observations of this type of deviation by CDF and DO [40) could be a hint that four 
generations with relatively light superpartners are present. 



8. Conclusions 

We have shown that a four-generation MSSM model is an attractive exten- 
sion of the usual three-generation MSSM in that gauge unification and automatic 
electroweak symmetry breaking via RGE renormalization both occur naturally for 
typical soft-SUSY-parameter boundary conditions at Mjj. A significant, but re- 
stricted range oft' and b' masses is allowed (see below) even after requiring that all 
Yukawa couplings remain in the perturbative domain throughout mz to Mjj evolu- 
tion; however, Yukawa unification is generally not possible. The domain of allowed 
t' and b' masses is impacted by constraints coming from the SUSY sparticle sec- 
tor, the most important of which is the requirement that the lighter r' eigenstate 
(t[) be more massive than the LEP-I limit of 45 GeV and also heavier than the 

LSP. This constraint can rule out a range of lighter b' masses; the smaller the 
soft-SUSY-breaking mo and rnj/2 parameters, the larger the range. For given t' 
and b' masses the r[ constraint places significant restrictions on mo and m^, and 

thence on other sparticle masses. Generally, small mo values are disfavored unless 
m~ is large, and there is a significant lower bound on m~ such that no mo, A choices 

yield a satisfactory r{ mass. There is no analogous bound in the three-generation 
case. Thus, the sparticle mass scale must be larger in four- vs. three- generation 
MSSM models. 



\ Note that Mg USY is not generally the same as the effective Msusy appropriate in discussing 
gauge coupling unification, which is sensitive to the entire particle content of the theory 
and not just the strongly interacting particles. 
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Significant differences are also found for the relations between different sparticle 
masses in comparing N g = 4 results to those for N g = 3. The most dramatic 
difference arises directly from the roughly factor of two larger value of the unified 
gauge coupling a(Mjj) in the N g = 4 case. This leads to ratios such as mj l m ^i 

and mj^/m~o being approximately a factor of two larger for N g — 4 compared to 
N g = 3 for a given choice of mo/m^. Indirect tests/verifications of the presence 
of a fourth generation by observing only a light gaugino and sparticles belonging 
to the first two families are thus possible. 

Direct experimental constraints on the perturbative four-generation MSSM 
model are becoming very strong, and near-term experiments could easily eliminate 
the model altogether. The v' and r' have masses that are strongly bounded from 
above for all wif, my > 45 GeV (the current LEP-I bound) and will be readily seen 
at LEP-176. For v' and r' masses just above the current LEP lower bounds (we 
take 50 GeV for these masses) perturbativity constraints force mp (pole) ^ 160 GeV 
and my (pole) < 120 GeV. We also can be quite certain that the t' is heavier 
than the b' since otherwise the t' decays to bW (or bW*), and t't' final states 
would have led to top-quark-like events at a rate inconsistent with CDF and DO 
results. For mp(pole) > my (pole), the b' can (and must, if it is to have evaded 
detection to date) decay via flavor changing neutral current processes to b + X 
(X = 7, g, qq, Z^*\ h°, . . .). For expected branching ratios, non-observation of b'b' 
production events at the Tevatron in which both b h s decay to b + 7 probably 
excludes my (pole) ^ 80—85 GeV. The range my (pole) > mz+my(pole)+ ~ 5 GeV 
(for which b' — > bZ decays have significant branching ratio) is probably excluded by 

the non-observation of the extra Z's expected from b'b' production events, except 
for soft-SUSY-breaking parameter scenarios such that m^o is small enough for 
b' — > bh° decays to dominate over b' — > bZ - - the scenarios for which the h° 
channel can be open have a large soft-SUSY-breaking scalar mass, mo. The top 
quark searches of CDF and DO would have detected production of t't' with t 1 —> b'W 
(or W*) and b' decaying hadronically, unless mp(pole) is quite close to my (pole) 
(implying soft W* decay products that evade their hard cuts). We roughly estimate 
that m t i(pole) must be ^ 115 — 120 GeV on this basis. All of this leads to a highly 
preferred mass range for the b' and t': 

120 GeV ;> m v (pole) Z my (pole) with 100 GeV ^ my (pole) ^ 80 GeV. (21) 

We estimate that, for an integrated luminosity of L = 150 pb~\ it is well within 
the capability of the Tevatron to either exclude all reasonable allowed models or 
detect a signal. 

Of course, some of the above restrictions follow from the assumption that the 
'top-quark' events at CDF and DO arise from the third- generation top. We have 
noted that the alternative case (emphasized in Ref. [14]) where these are t't' events 
in which t' — * bW, while the normal t has mt ~ mm and decays via the difficult- 
to-detect t — > t\Xi channel is not possible in the SUSY GUT context assuming 
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universal mo and A parameters. As one adjusts parameters so as to minimize 
m^r + m~o, the t[ always becomes too light, its mass falling below the LEP-I 

bound or below the LSP mass, or the theory becomes inconsistent by virtue of the 
CP-odd Higgs or some colored and/or charged 4th-family sparticle being required 

to have m 2 < 0, long before the lighter t\ reaches masses below rayy. 

Certainly, if there is a fourth generation with Yukawa couplings that remain 
perturbative up to Mjj, experimentalists will discover a plethora of new signals at 
LEP-II and with increased luminosity at the Tevatron. If in the end no signal is 
found, requiring (roughly) m T >, m v i ^ 95 GeV and , my ~ 200 GeV, determined 
four-generation model builders must become resigned to having one or more of 
the Yukawa coupling constants becoming non-perturbative before evolution up to 
Mjj is complete. While there is no known fundamental reason to disallow this, 
such a scenario is distinctly less predictive, and therefore less attractive, than full 
perturbative evolution up to the unification scale for all parameters. 

Acknowledgments 

This work has been supported in part by Department of Energy grants #DE- 
FG03-91ER40674, #DE-FG02-85ER40214, and by the Davis Institute for High 
Energy Physics. 

APPENDIX 



In this Appendix we collect all of the renormalization group equations for 
the four-generation extension of the MSSM that are used in our calculations. 
Among the many references available, we primarily employed Bjorkman and Jones, 

Ref. [15], Martin and Vaughn) 411 and Cvetic and Preitschopf 161 for the various beta- 
functions required. 

Assuming no KM mixing and universal scalar masses at Mj/, and defining t = 
(^)log[Q( GeV)], yi = Xf/(4ir), and a{ = gf/(4ir), we summarize the RGE's 

below. We denote the anomalous-dimension contributions to by An'y^ at one 

loop and (47r) 2 7^ at two loops for the i th (unmixed) fermion. 

For the Yukawa couplings we write 

d f t = + (4^7? } ) , (A.l) 

where i = t,b, r, t', b', is', r' . The expressions for A.ii^p and (Air) 2 ^ 2 ^ are as follows: 
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(l) 13 16 
4^ = ~J^ a i ~ 3 «2 - y«3 + vyt + 3j/f +yb + Vv 

,„ x2 (2) 3523 o 27 2 80 2 136 
(4tt) 7t v = yy«i + y«2 + y«3 + «i«2 + "^-«i«3 

6 4 

+ 8a 2 «3 + 1%«3 + 1 6^0:3 + + -yt>ai 

+ \vbOt\ + Qyt®2 - 22y 2 - 9y\, - 3y 2 v , - 5y b y t 
5 

- Zyvyv - y^yv - 5y 2 - 9ytyt> - tybyv - yby^ 

- yby T ' - ?>ytyv' ; 



4iry b = -yjr«l - 3a 2 - y «3 + %& + 3y & > + yt + y r + y T 

, 2 ( 2 ) 371 2 27 2 80 2 8 
■x) % = -gQ-ai + y«2 + y"3 + «i«2 + g«i«3 

2 4 

+ 8a 2 «3 + l%b«3 + 16yfe'a 3 + -y h a\ + -yt®i 
2 6 6 

- 72/6' a l + 7^1 + + %&«2 
5 5 5 

- 22^ - 9y 2 , - 3y 2 - 3y 2 T , - 5y b y t - 3y b ,y v ~ y^u' 

- 5yf - 9y b ?/ & / - 3y t y t > - 3y b y r - 3y 6 y T / - y t y v ' ; 



47T7r = --ai - 3a 2 + 3y 6 + 3y & / + 4y r + y T > 
5 

(4tt) 2 7| 2) = —a? + — a 2 + T«ia2 + 16y6«3 + 16yb>a3 
10 z 5 

2 2 6 6 „ 

- 72/6ai - ^'"l + 7Vr a l + -y T '<X\ + %t"2 
5 5 5 5 

- 9y\ - 9yl - 10y 2 - 3y 2 , - 3y b y t - Zy v y v - y T >y u ' 

- 3y T y T ' - 9y T yb - Vyryv ; 
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,„ n2 (2) 3523 o 27 2 80 2 136 

(47r) 7., = cti H a 9 H Qo + mao H 0110:3 

v ) H 450 1 2 2 9 3 45 

6 4 

+ 8a 2 a 3 + 16^a 3 + I6y t as + -yt'Oii + -ytai 



+ ^y&'Oi + 6^o 2 - 22y 2 - 9y\ - 3y 2 , - hy b >y v 


- 3y & y 4 - y T 'y„' - 5y 2 , - 9y t >y t - 3y b ,y b - y b ,y T > - y b ,y T - 3y v y v > ; 

(A.5) 



7 16 

4^7^ = - 3«2 - y«3 + fyb* + 3y b + y v + y T + y T > 



,2 (2) 371 2 27 2 80 2 8 
(4tt) 7^ = — a x + y« 2 + — o 3 + oio 2 + -aia 3 

2 4 

+ 8a 2 «3 + 16yfe'a 3 + 16y b o 3 + -y b >ai + -y*'Oi (^.6) 

2 6 6 „ 

- rVbUi + rVr' a l + tVt^i + 6y b >a2 

5 5 5 

- 22yl - 9y% - 3y 2 T , - 3y 2 - 5y v yt> - 3y b y t 

- y T 'yu' - 5y 2 ' ~ $VMb - ZytW - 3y b >y T - 3y b <y T > - y v y V ' ; 



d) 9 

47T7;/ = --01 - 3a 2 + 3y v + 3y b + y T + Ay T > + y v , 

X 71 27 9 

(4tt) 2 7| 2) = - jQ-ai + y«2 + ^ia 2 + 16y 6 /o 3 + 16y fe o 3 

2 2 6 6 „ , . 

- -y b 'Oi\ - -y b a.\ + -y T 'a\ + -y T ai + 6y r 'a2 (A.7) 
5 5 5 5 

- 9y 2 / - 9y 2 - 10y 2 , - 3y 2 - 3y v y t > - 3y b y t 

- y T <yv> - 3y T y T > - 9y T >y b > - 9yr>y b - 3y v > 

- 3y v *yt - 3„>y t > - 2y T >y v > ; 
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(l) 3 

47 H> = - 3a 2 + i>yv + 3yt + 4jv + y T > 

, A n2 (2) 267 2 27 2 9 ,„ ,„ 6 

(4tt) 7^/ = -gQ-ai + y«2 + gtti«2 + 16y t a 3 + 16?/*' a 3 + -y r '«i + 6tv« 2 

6 4 4 2 2 2 

+ -=yv>OL\ + -y t ai + -yt>a\ - 9y t - 9y t > - Wy^ ~ 3y b y t - ZyvW 



- tyr>vv - Zyl> - y T yr' - Zy T >yb - ZyT>yv - 9y v >y t 



For the Higgs mass terms, we have: 



(A.8) 



dTf? 3 

= 3 J2 y^i + y» ,u "' - - 3 «2^ 2 ; (a.9) 



dt 

q=t,t' 



q=b,V 1=t,t' 

The third and fourth generation scalar quark and lepton soft mass terms evolve 
according to: 

= yt,t>U t ,t> + y b ,b>Vb,b> ~ ^mMf - 3a 2 M| - y« 3 M 3 2 ; (A.ll) 



dlTtr r tji 



U,U> _ 7 , 16 2 16 71^2 



(A.12) 
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(A.13) 



dm? r, 



i > 3 

^— = yr,T'V r ^ + b iL ,y v ,U v i - -a\M\ - 3a 2 M% ; 



(A.14) 



rfm l E> _ 12 2 

— = 2y T;T ,V T!T , - — ai M?; = 2yM< . 



(A.15) 
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In the above, 



u t,t> = m 2 H 2 + m Q,Q> + m u,v + A lt>, 



(A.16) 



U v * = m 2 H2 + m 2 L , + m 2 N + A 2 ,, (A.17) 
= m 2 Hl + m| Q , + D , + (A.18) 

The running of the first and second generation scalar quark and lepton soft-SUSY- 
breaking masses is obtained by neglecting all Yukawa couplings in the expressions 
above. The parameters fx, specifying mixing of the Higgs superfields in the super- 
potential, and B, such that Bfi is the coefficient of the H1H2 scalar field mixing 
term in the soft-SUSY-breaking potential, evolve according to: 



dlog(fi) 3 



I y i + \ ^-7^1-^2; (A.20) 



dt 

q=t,t',b,b' 1=t,t',V 

dB 3 

— = 3 Yl yi A <i + Yl yi A i + -a 1 M 1 + 3a 2 M 2 . (A.21) 

q=t,t',b,b' l=T,r'y 

The running of the third and fourth family soft-SUSY-breaking potential tri-linear 
term A coefficients is given by: 



=6A t:t >yt,t> + A bjb ,y bjb , + 3A t > :t yt',t + A v ty v , 

13 16 
+ — aiMi + 3a 2 M 2 + —a 3 M 3 ; 
15 3 



(A.22) 



dA b , 



'- =G\b'Vb,b' + M,Wt,v + 3A h ^ b y b , )b + Y A iVi 



dt 

i=TS (A.23) 

7 , r „ , , 16 
+ —aiMi + 3a 2 M 2 + — a 3 M 3 ; 
15 3 

dAq- ft ^ — "\ 
dt =4:A T yy T y + A T ^ T y T ^ T + 3 2^ AqVq + k,T>A v iy v i 

q=b,b' (A. 24) 

9 

+ - ai Mi + 3a 2 M 2 ; 
5 
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—± = AA v ,y„, + A T ,y T , + 3 ^ A q y q + 3a 2 M 2 + - ai Mi . (A.25) 

«=*,«' 

The Higgs boson, scalar quark and scalar lepton mass-squared RGE's generally 
include a term in each equation that is a numerical factor times a\S' , where 



(A.26) 



S' = m% 2 - m 2 Hl +rriQ + m 2 QI - 2m^ + m 2 D 
— 2m 2 j l + m 2 D , — m 2 L — m 2 L , + m 2 E + m 2 E , . 

Substituting from the RGE's with S' terms included, one finds 

W = T aiS ■ (A ' 27) 

Because of the universal scalar mass boundary condition at Mjj that we assume, 
S'(Mjj) = 0. Therefore S'(t) = for all t and we do not include this term in the 
mass-squared RGE's presented above. 

The final soft-SUSY-breaking parameters are the gaugino masses. Using the 
notation N g for the number of generations, these evolve according to: 

dM 2 



dt 
dM s 



= (2N g -5)a 2 M 2 ; (A.28) 



(2.V,, !))o ;! .U ; .. 



Finally, we present the gauge coupling RGE's for completeness. These are the 
same as the massless fourth generation neutrino case, since the extra right handed 
neutrino that provides a mass term in our case is decoupled from the gauge sector. 
Keeping the number of generations, N g , explicit, including the third and fourth 
(unmixed) generation Yukawa couplings, and including exactly two Higgs doublets 
(as appropriate for the MSSM), one has the gauge RGE's: 



dai 
~dT 



= (2N g + °-}a 2 + ^ 
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38 ir 9 

— N a H 

15 9 25 



18, , 14. . 26, 

-^{Vr + Vt') - —{yb + Vb') - y {yt + y? 



6 Ar 9\ 88 Ar 
-N g + - )a 2 + -N g a 3 



(A.29) 
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(^N g + ^ «i + (14A^ 5 - 17)a 2 + 8N s a 3 



- 2(y r + y T /) - 6(y b + y & , + y f + yt') 



(A.30) 



da% 
~df 



(2N g - 9)a 2 3 + ^ 



N g ai + 3N g a 2 + (^j N g ~ 54^ 

- Kvb + yy + yt + w 



11 

15" 



(A.31) 



The renormalization group equations were implemented numerically, and iter- 
ation was employed to find a fully consistent solution for the complete evolution 
between Mjj and mz- 



APPENDIX 

In this appendix we discuss the question of the accuracy with which cos 2/3 
and, thence, tan/3, can be determined using the mass differences of Eq. (10). 1421 

To use m~ —m~ requires accurate determinations of m~ and tuj^. The required 

measurements are best performed at an e + e~ collider where the center-of-mass 
energy is precisely known and the energy spectra of final leptons can then be 
directly related to the masses of the produced particles and their decay products. 

Determination of mj^ is relatively straightforward since £l — > is almost always 
an important, if not dominant, mode — if mj > m~+ ~ m x°> then — > Xi u 
and £l — > ^£ can be competitive. A fit to the i spectrum will yield the masses 

of both the £l and the Xi- Should the X\ v an( i modes also be important, 
the fit would presumably reveal the presence of several upper end points, and 
allow determination of both m~o and m~o. A Monte Carlo study is desirable to 

Al A2 

determine the exact accuracy with which all masses could be determined in the 
multiple decay case. Perhaps accuracies of about ±10 GeV could be achieved. For 
mj ~ 240 GeV (near our lower limit in the N g = 4 case scenarios, but already 

rather marginal for a y/s = 500 GeV NLC) this would correspond to accuracies of 
~ 4%, i.e. on the edge of what would be useful. 

Determination of m~ follows a similar pattern. We give below discussions for 
two scenarios: a) m~ < m~+ and b) m~ > m~+ . Scenario a) does not arise for the 
specific N g = 4 models illustrated in Figs. 13-15, but might be relevant if there 
was significant non-universality in the Mjy-scale soft-SUSY-breaking parameters. 
For N g = 4, universal ttiq and Mj values always lead to m~ > m-p as shown. 
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If m~ < m~+ then the xt decays to a mixture of £^v and lv\ since m~ > m~ 

Al ^ &L 

the latter decay is usually the dominant one — indeed, in many cases mj — m~ 

is large enough that the xt — *■ £l v channel is forbidden even though the xt ~^ & 
channel is open. At an e + e~ collider the end-points of the £ spectrum from xtxi 
production followed by xt & decay provide a determination of both m~+ and 

Al 

m~; direct observation of the v is not required, a fortunate fact given that it decays 
invisibly to vx\- 

If m~ > rn~+, then v — > X\^ an d the xt decays via a virtual sneutrino or 

slepton to the three-body mode lvx\, via virtual squark to qqx\i or via virtual or 
real W to a mixture of these two final states. Analogous to the previous case, the 
£ spectrum end-points will allow determination of m~ and m~+ , although one will 

have to carefully account for the underlying smooth £ spectrum from the xt decay 
if events containing only 2£ + 4j+ missing energy (which avoid contamination from 
the soft £(s) from xt decay) do not occur at a sufficient rate. Ideally, it would be 
good to learn first about the xtxi channel by setting y/s to a value above 2m~± 
but below (or not far above) 2m~ — due to the rapid turn-on of a two-fermion 
channel, the rate for xtxi production might be reasonable, while the uu rate could 
be zero or quite small (due to a large — m~+ splitting and/or the slow turn-on 

of the vv spin-0 pair channel). To study vv, one would then up the energy to a 
level such that the vv rate was large. 

A significant complication is the fact that £l^l production would also be 
present at a similar rate to vv (recall that mj^ m~). For the present case of 

m~ > rn~±, £l — > £xi, u Xi^X2 wm a ^ occur > most probably with similar branch- 
ing ratios. Events in which both £^s decay to £xi can De eliminated by requiring 
some jet activity (which is generally present for the vv final states). If both i^s 
decay to X\ v i then the only charged Ps would be from decays of the xf's, which 
would yield a smooth (and relatively soft) £ spectrum; in addition, these £l£l 
events would tend to have larger missing energy than the events. If one £l 
decays to £x\ and the other to vxt \ followed by xt — >• £vx\i then we could end up 
with a 2!+ 2j+missing energy final state that might not be easily distinguished 
from the 2£ + 4j+missing energy state of interest for vv events (certainly not all 
jets in the latter case would be detectable). But, the spectrum of the £ from the 
xt decay would be soft and smooth, so that the ^-spectrum thresholds present in 
^ ~" *■ ^X.t m ight still stand out. 

Thus, it would seem that there are means for isolating the vv events of interest 
and extracting the £ spectrum thresholds that would allow a reasonable deter- 
mination of m~ (especially if m~± is known from a lower energy measurement). 

Al 

Unfortunately, these complicated scenarios have not been explicitly studied for 
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an e + e~ collider. The results for simpler scenarios (see Ref. [30] and references 
therein) suggest the crude estimate that mass determinations for m~ might be 
possible within ±5%, at least for masses of order 200 GeV and below {i.e. compa- 
rable to the minimum possible m~ values in the scenarios discussed in the previous 
sub-section.) However, including the fact that the mj and m~ errors must be 

combined in quadrature, despite this fairly small uncertainty in the mass measure- 
ments we would not have the accuracy required for a 3<r sensitivity to cos 2(3, even 
for the smallest possible masses allowed in the N g = 4 case. 

Let us now turn to the m~ — m~ measurement. First, there is a very real 

possibility (a near certainty for N g = 4) that the squarks are simply so heavy, 
m~ ^ 300 GeV, that impossible accuracy would be required for sensitivity to cos 2(3. 
Even for masses below this level, it is still far from clear that the required accuracy 
can be achieved. Squarks will decay to a quark plus real or virtual gluino, the latter 
then decaying to a variety of final states (with ggx? unlikely to be dominant for the 
larger m~ predicted for N g = 4). Jet spectra end-points (which would appear on 
top of a smooth jet-spectrum background from the real or virtual g decays) could 
provide a certain level of accuracy in the squark mass determination, but there 
would be a lot of overlap of the spectra from different g decay channels, and of the 

spectra from the d and u decays themselves, that might very well prevent a good 
determination of m^—m~. An extremely careful Monte Carlo study is required to 
be able to determine with any certainty the level of accuracy that can be achieved. 

Overall, we conclude that our ability to determine cos 2(3 and thence tan (3 
must remain a topic of further study. Fortunately, the uncertainty in cos 2(3 is not 
the limiting factor determining the level of accuracy with which other mass sum 
rules and relations can be tested. Other experimental and theoretical uncertainties 
are almost certainly more important. 



•52 



REFERENCES 



1. J.F. Gunion, D.W. McKay and H. Pois, Phys. Lett. B334, (1994) 339. 

2. For a detailed discussion of the neutrino sector, see S. King, Phys. Lett. 
B281, 295 (1992). 

3. OPAL Collaboration, M.Z. Akrawy et al, Phys. Lett. B236, 364 (1990); 
ALEPH Collaboration, D. Decamp et al, ibid., 236, 511 (1990). A recent 
survey of lower bounds on heavy and neutral charged leptons and on heavy 
quarks is given by the Particle Data Group, K. Hikasa et al, Phys. Rev. 
D50, Part 1 (1994). 

4. M. Drees and S.P. Martin, 'Implications of SUSY Model Building', Report of 
Subgroup 2 of the DPF Working Group on "Electroweak Symmetry Breaking 
and Beyond the Standard Model" , MADPH-95-879. 

5. H. Baer, J. Gunion, C. Kao and H. Pois, Phys. Rev. D51, 2159 (1995). 

6. An exception is the early N g = 4, m v > = supergravity extension of the SM 
by M. Cvetic and C. Preitschopf, Nucl. Phys. B272, 490 (1986). 

7. For more recent work, see P. Agrawal and W. -S. Hou, Phys. Rev. D46, 
1022 (1992); P. Agrawal, S. Ellis and W. -S. Hou, Phys. Lett. B256, 289 

(1991) ; S. King, Ref. [2]; C. Hill and E. Paschos, ibid., 241, 96 (1990); 
W.-S. Hou and R. G. Stuart, Nucl. Phys. B349, 91 (1991); M. Sher and 
Y. Yuan, Phys. Lett. B285, 336 (1992); H. Fritzsch, Phys. Lett. B289, 92 

(1992) ; B. Mukhopadhyaya and D.P Roy, Phys. Rev. D48, 2105 (1993); J.L. 
Hewett, SLAC-P UB-6521 (1994); Z. Berezhiani and E. Nardi, UM-TH-94-36, 
|hep-ph/9411249| (1994). 

8. Mark II Collaboration, G. S. Abrams et al, Phys. Rev. Lett. 63, 2173 
(1989); L3 Collaboration, B. Advera et al, Phys. Lett. B231, 509 (1989); 
OPAL Collaboration, I. Decamp et al, ibid., 231, 519 (1989); DELPHI 
Collaboration, M. Z. Akrawy et al, ibid., 231, 539 (1989). 

9. For a recent example of theoretical motivation of N g = 4, see Berezhiani and 
Nardi, Ref. [7]. 

10. P. Langacker, in Proceedings of the PASCOS-90 Symposium, eds. P. Nath 
and S. Reucroft (World Scientific, 1990); P. Langacker and M. -X. Luo, Phys. 
Rev. D44, 817 (1991); U. Amaldi, W. de Boer, and H. Furstenau, Phys. 
Lett. B260, 447 (1991); J. Ellis, S. Kelley and D. V. Nanopoulos, ibid. 260, 
131 (1991). 

11. P. Langacker and N. Polonsky, Phys. Rev. D47, 4028 (1993); ibid., 49, 1454 
(1994); V. Barger, M. Berger, P. Ohmann, and R. J. N. Phillips, Phys. Lett. 
B314, 351 (1993); G. Kane, C. Kolda, L. Roszkowski, and J. Wells, Phys. 
Rev. D49, 6173 (1994); G. Ross and R.G. Roberts, Nucl. Phys. B377, 571 
(1992); B.D. Wright, MAD/PH/812. For string Yukawa threshold effects, 
see I. Antoniadis, E. Gava, K.S. Narain and T.R. Taylor, Nucl. Phys. B407, 
706 (1993). 



53 



12. P.H. Frampton, S. Nandi and J.G. Scanio, Phys. Lett. B85, 225 (1979). 

13. S. Dimopoulos and A. Pomarol, CERN-TH/95-44, |hep-ph/9502397| . 

14. H.E. Haber, preprint CERN-TH/95-178 (1995). 

15. D. V. Nanopoulos and D. Ross, Phys. Lett. B118, 99 (1982); J. Bjorkman 
and D. Jones, Nucl. Phys. B259, 533 (1985). 

16. P. Langacker and N. Polonsky, UPR-0642T 

17. M. L. Swartz, SLAC Report No. PUB-6710 

18. J. Bagger, V. Matchev and D. Pierce, Phys. Lett. B348 443 (1995). 

19. See, for example, V. Barger, M. S. Berger and P. Ohmann, Phys. Rev. D47, 
1093 (1993). 

20. F. Abe et al. (CDF Collaboration), Phys. Rev. Lett. 74, 2626 (1995); S. 
Abachi et al. (DO Collaboration), Phys. Rev. Lett. 74, 2632 (1995. 

21. B. Pendelton and G. Ross, Phys. Lett. B98, 291 (1981); C. Hill, Phys. Rev. 
D24, 691 (1981); W. Bardeen, M. Carena, S. Pokorski and C. Wagner, Phys. 
Lett. B320, 110 (1994); J. Bagger, S. Dimopoulos and E. Masso, Phys. Rev. 
Lett. 55, 920 (1985). 

22. See, for example, G. Gamberini, G. Ridolfi and F. Zwirner, Nucl. Phys. 
B331 331 (1990). 

23. T.L. Barklow et al., 'Electroweak Symmetry Breaking and Beyond the 
Standard Model', SLAC-PUB-95-6893 (1995). 

24. J.F. Gunion, A. Stange and S. Willenbrock, in preparation. 

25. S. Mrenna and G. Kane, CALT-68-1938 (1994). 

26. These effects were first explored in J.F. Gunion and S. Geer, Proceedings of 
the "Workshop on Physics at Current Accelerators and the Supercollider'' ', 
eds. J. Hewett, A. White, and D. Zeppenfeld, Argonne National Laboratory, 
2-5 June (1993), ANL-HEP-CP-93-92, p. 335. 

27. M. Drees and M.M. Nojiri, Phys. Rev. D45, 2482 (1992). 

28. S. Martin and P. Ramond, Phys. Rev. D48, 5365 (1993), provide a convenient 
summary of squark and slepton mass relations in the N g = 3 case. 

29. J.F. Gunion and H.E. Haber, Phys. Rev. D37, 2515 (1988). 

30. See H. Baer et al., "Low Energy Supersymmetry Phenomenology", preprint 
FSU-HEP-950401, and references therein. 

31. V. Barger, M. Berger, J. Gunion and T. Han, preprint UCD-95-12 (1995). 

32. F. Abe et al, Phys. Rev. Lett. 64, 147 (1990), Phys. Rev. Lett. 68, 447 
(1992); T. Trippe, Particle Data Group, private communication. 

33. J. Houston and R. Blair, private communication. 



54 



34. W.-S. Hou and R. Stuart, Phys. Rev. D43, 3669 (1991). 

35. D. Schaile, talk presented at the 27th Int. Conf. on High Energy Physics,, 
Glasgow, Scotland, July 20-27, 1994; SLD Collaboration, F. Abe et al, Phys. 
Rev. Lett. 73, 25 (1994). 

36. M. Peskin and T. Takeuchi, Phys. Rev. D46, 381 (1992). 

37. C. P. Burgess, S. Godfrey, H. Konig, D. London and I. Maksymyk, Phys. 
Rev. D49, 6115 (1994); C. P. Burgess, S. Godfrey, H. Konig, D. London and 
I. Maksymyk, Phys. Lett. B326, 276 (1994); I. Maksymyk, C. P. Burgess and 
D. London, Phys. Rev. D50, 529 (1994); P. Bamert, C.P.Burgess MCGILL- 
94-27, (1994). 

38. A. Kundu and P. Roy, SINP-TNP/94-07, TIFR/TH/94-38. 

39. H. Pois, T.C. Yuan, unpublished 

40. Figures available on the CDF WWW Home Page and private communica- 
tions. 

41. S. Martin and M. Vaughn, Phys. Rev. D50, 2282 (1994). 

42. We thank H. Murayama for several conversations on various aspects of the 
discussion in this Appendix. See also Ref. 30. 

FIGURE CAPTIONS 

1) Gauge coupling unification, including two-loop gauge contributions but not 
two-loop Yukawa contributions; dashed (solid) curves correspond to N g = 
3 (4). The illustration is for sin 2 6 W = 0.2316 and a^(m z ) = 127.9. 

2) We give the allowed rrif, my parameter space regions for m t (m t ) = 165 GeV 
in the cases tan (3 = 1.5 and 2.2. We have taken m v > = m T > = 50 GeV. Full 
two-loop contributions to the gauge coupling beta functions are included. 
Small dots indicate regions of ,my disallowed by demanding perturbative 
Yukawas (A^ < 3.3) and gauge coupling unification to within 2.5%. Small 
squares indicate additional points excluded if gauge unification is required to 
better than 0.01%. 

3) Contours in my and mp (running mass) parameter space of constant Mjj and 
otz(mz) after the inclusion of Yukawa terms in the two-loop gauge coupling 
beta functions, assuming mt(mt) = 165 GeV and tan/3 = 1.5. We have taken 
m v i = m T i = 50 GeV. 

4) We display contours of constant Xy, Xf and X T > at Mjj in the region of mp, my 
parameter choices allowed for 2.5% unification accuracy with m t {m t ) = 
165 GeV and tan/3 = 1.5 or 2.2. We have taken m v > = m T i = 50 GeV. 
Two-loop contributions to the gauge coupling beta functions are included. 

5) In a) we illustrate the evolution of the soft-SUSY-breaking scalar mass- 
squared parameters m 2 H and m 2 ^ as well as the value of In b) we show 
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the evolution of the coefficients mf, m\ and ml of the Hf, iff and C 1E 1 E 2 
scalar field potential terms, as well as that of S = m\m\ — m\. We have 
taken £ = l/v 7 ^, U = -1, m 1/2 = 400 GeV, tan/3 = 1.5, = 100 GeV, 
m^/ = 100 GeV, and m T > = m v > = 50 GeV. All masses are given in units of 

6) Contours of /x, m^o and m^o in the me, parameter plane for tan/5 = 1.5, 
m t (mt) = 165 GeV and m^ 2 = 600 GeV. We give results for the Mjj-scale 

dilaton boundary conditions: £o = l/y/3, — ~ 1- The diamonds indicate 

regions of parameter space disallowed because m~/ < max |45 GeV, j at 

energy scale m^. 

7) We plot the ratio of N g = 4 to N g = 3 values for: F(h -> gg),T(h -> 77), 
and r(/i — > gg) x BR(h — > 77). We adopt the scenario of tan/? = 1.5, 
mj = 165 GeV, = = 100 GeV, r?v = 77v = 50 GeV, and assume 
superpartners are sufficiently heavy that their contributions to these one-loop 
quantities are small. 

8) Evolution of m~, , m~ , and m~, (in units of m? /0 ). We have chosen tan 3 = 

1.5, £ = l/v 7 ^, £a = -1, and m 1/2 = 600 GeV. Results for m v = 100 GeV 
and my = 90 GeV are compared to those for m t < = 100 GeV and my = 
110 GeV. The corresponding m^-scale m~^ values are 189 and 311 GeV, 

respectively. 

9) We display the r{ constraint boundaries for a variety of Mjj-scale scenarios. 
Regions to the left of the boundaries are disallowed. We have taken m% = 
165 GeV, 77v = rn V ' = 50 GeV, and tan 8 = 1.5. 

10) We plot the minimum value of £0 that is allowed by the constraint m~i > 

max 1 45 GeV, m~o j for a given value of m^ 2 after scanning over all possible 

values of A, assuming running masses mt = 165 GeV, m T > = m v > = 50 GeV, 
my = m? = 100 GeV, and tan (3 = 1.5. 

11) We plot the minimum possible value of m~(m~), and the corresponding values 
of m^o and my (pole) —m^o — m^, as a function of £4 for a series of £0 values: 
£o = (solid); £0 — 0.5 (long dashes); £0 = 1 (dots); £0 — 2 (dot-dash); 
£o = 3 (short dashes); £0 = 5 (dash-dot-dot). We have taken m t (m t ) = 
165 GeV, my(my) = m t >{mt>) = 100 GeV (corresponding to pole masses of 
approximately 175 GeV and 105 GeV). The minimum m~ is that allowed by 
the t[ constraint for a given £o>d choice. The corresponding tan/3 value is 
either 1.5 or 1.6 in all cases. At each £0 an d £a value, all values of m^ 2 < 
10 TeV were scanned. Curves terminate when no consistent solution is found. 

12) Mass spectra for the t[ for the three standard (£o,d) choices. At the lower 
values of m~ the spectrum terminates either because m~, falls below 45 GeV 
(high-mo scenario), or because m~i < m~o (dilaton and no-scale scenarios). 
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We have scanned in tan/3, fixing mt(mt) = 165 GeV, m v > = m r > = 50 GeV 
and m? = my = 100 GeV. 

13) Mass spectra for dilaton boundary conditions. Masses are given in units of 
m~. We have taken m T > = m v > = 50 GeV, my = my = 100 GeV, /i > 0, and 
scanned in tan /3. 

14) Mass spectra for no-scale boundary conditions. 

15) Mass spectra for the large mo boundary condition case. 

16) A plot of r~ and rj (see text) as a function of £o for the typical case of 
m~o = 117 GeV. We have taken tan/3 = 1.5, mt = 165 GeV, and, in the 

Al 

N g = 4 case, m? = my = 100 GeV, m T > = m v > = 50 GeV. 

17) We plot the maximum value of V^h for which T(t' — > b'W*) > F(t' — > bW) as 
a function of m? for fixed values of my = 80 GeV and 100 GeV. All masses 
are pole masses. 

18) We show the 95% C.L. and 90% C.L. constraints on the my (my) , m t > (mf) 
parameter space (assuming m T > = m v * = 50 GeV and tan/3 = 1.5) from 
Z-pole precision electroweak measurements (coupled with mw/mz). We 
consider two reference SM's specified by mt(m t ) = 165 GeV, mjj = 100 GeV, 
and the two values, a s (mz) = 0.12 (dashed lines) and 0.13 (solid lines), see 
text. 

19) We plot the percentage enhancement R as (see text) as a function of Q for 
N g = 4 choices of = 200 GeV and M$°sy = 400 GeV - 
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Excluded Regions: m t (m t ) = 165 GeV, rn b (m b )=4.6 GeV 
Unification to within 2.5% (.01%): •■•(••.+ □□) 
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Fig. 3 
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m b =4.6, m t =165, m t . = m b .= 100, m„, = m T , = 50 GeV 
tanj8=1.5, f =l/V3, f A = -l, m 1/2 = 600 GeV 
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m b = 4.6, m t =165, m„. = m T . = 50 GeV 
tan/S = 1.5, f =l/V3, £ A =-1, m 1/2 = 600 GeV 




80 90 100 110 80 90 100 110 120 



m b . [GeV] m b . [GeV] 



Fig. 6 

59 



Hiess Detection at the LHC: N =4 vs. 3 
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m b = 4.6, m t =165, m t .= 100, m l/ . = m T , = 50 GeV 
tan/S=1.5, f =l/V3, f A =-l, m 1/2 = 600 GeV 
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Minimum £ for Allowed Solution 
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Results at Minimum Possible m 
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No-Scale Mass Spectra: N =4 
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High — m Mass Spectra: N — 4 
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m b = 4,6, m t =165, m t . = m b .= 100, m 1/ . = m T . = 50 GeV 
tan/S=1.5, f A = 0, m 1/2 = 600 GeV 
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Excluded by Precision EW at 90% and 95% CL 
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